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In  section  2,  we  calculate  surface  wave  magnitudes  for  1 12  Nevada  Test  Site  (NTS)  explosions  from  a  data  set  of 
regional  long-period  seismograms  from  North  American  stations.  In  order  to  utilize  the  nearer  regional  stations  (^<25*^),  a 
new  method  for  determining  Ms  and  been  developed  which  employs  synthetic  seismograms  to  establish  a  relationship 
between  the  amplitude  of  the  regional  Airy  phase,  or  Rayleigh  pulse,  of  the  data  and  an  associated  surface  wave  magnitude, 
based  on  conventional  Ms  determinations,  calculated  form  a  synthetic  seismogram  propagated  to  40*’.  This  method  easily 
lends  itself  to  implementing  path  corrections.  The  inclusion  of  path  cotrcctions  decreases  the  Mj  variance  by  a  factor  of 
two  and  also  decreases  the  average  value  by  0.08  magnitude  units.  This  latter  effect  is  attributed  to  the  particular  station 
network  used.  The  method  gives  stable  M^  values  that  correlate  well  with  other  magnitude  scale  values  over  a  range  of  three 
orders  of  magnitude  in  source  yield.  Our  most  refined  M^  values  yield  the  relationship  mb  =  0.9  x  Mg  t  1 .82.  When  events 
are  grouped  with  respect  to  source  region,  significantly  better  fits  to  these  individual  site  linear  regression  curves  are 
obtained  compared  to  the  fits  obtained  using  a  single,  all  inclusive  model.  This  observation  implies  that  shot  site 
parameters  and  source  structure  effect  surface  wave  magnitude  measurements,  although  event  yield  site  distribution  also  may 
be  in  part  responsible. 

Since  our  magnitude  values  arc  based  on  a  theoretical  continental  structure,  we  regressed  our  values  with  the  more 
standard  values  of  Marshall  et  al.  (1979).  Using  14  common  NTS  events  we  found  that  our  values  were  greater  by  0.53  +  - 
0.03  magnitude  units. 

In  section  3,  we  present  a  note  on  the  relocation  of  Tibetan  earthquakes.  Present  ISC  locations  for  earthquakes 
beneath  Tibet  indicate  a  random  distribution  of  events  down  to  a  depth  of  about  SO  km.  This  distribution  might  be  expected 
from  a  relatively  cool  crust  which  would  allow  the  scismo-gcnic  zone  to  extend  to  such  depths.  A  detailed  investigation  of 
the  Tibetan  earthquakes,  with  magnitude  greater  than  5.5  from  1964  to  1986,  yields  a  distinctly  different  picture. 

Waveform  modeling  of  depth  phases  such  as  pP  indicates  that  only  three  or  four  events  from  this  population  is  actually 
deeper  than  25  km.  These  few  events  occur  near  the  edges  of  the  Plateau  where  active  subduction  is  occurring  as  suggested 
by  the  thrust-like  nature  of  their  mechanisms.  The  events,  averaging  the  entire  population,  occurred  earlier  than  indicated 
by  the  (SC  by  about  3  seconds  which  leads  to  about  a  1.5%  and  0.5%  over  estimation  of  and  Sn  velocities  respectively 
applying  ISC  tables  and  standard  flat-layered  models.  A  more  serious  error  occurs  if  the  P^  and  Sn  velocities  are  determined 
by  correcting  for  source  depth  but  assuming  the  ISC  origin  times. 
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Summary 


The  research  performed  under  the  contract,  during  the  period  1  August  1989 
through  31  July  1991,  can  be  divided  into  three  main  topics;  modeling  regional 
broad-band  seismograms  from  the  Imperial  Valley  to  Pasadena,  determining  surface 
wave  magnitudes  for  NTS  events  using  regional  data,  and  waveform  modeling  to 
determine  the  depth  of  earthquakes  in  Tibet  leading  to  new  origin  times  and  a  better 
estimation  of  and  Sn  velocities. 

In  section  1,  we  address  broad-band  wave  propagation  along  a  corridor 
from  Imperial  Valley  to  Pasadena,  California.  The  path  consists  of  50  km  of  slow 
basin  structure  with  a  shallow  moho  followed  by  250  km  of  relatively  normal 
Southern  California  structure  approaching  Pa.sadena.  Rvents  occurring  in  the  valley 
produce  extended  body  wave  codas  which  carry  a  disixirsive  wavefonn  imprint 
along  with  longer  period  Rayleigh  waves.  Many  features  of  these  seismograms  at 
period  greater  than  a  few  seconds  can  be  model^  by  applying  a  finite-difference 
technique  to  a  2D  structure.  Shorter  periods  display  complex  behaviors  not  easily 
model^  by  present  techniques,  especially  for  shallow  events.  Events  occurring  in 
the  normal  structure  produce  broad-band  SH  waveforms  that  can  be  modeled, 
analytically,  if  we  include  some  extra  empirically  derived  parameters  determined  by 
calibration  events.  TTie.se  parameters  are  determined  by  adjusting  travel  time 
differentials  between  the  decomposed  wavefield  S,  sS,  etc.  For  example,  delays  in 
sS  relative  to  S  correspiond  to  local  corrections  for  sediment  cover  for  a  particular 
source  region.  Asperities  and  directivity  become  increasingly  important  at  shorter 
periods,  as  demonstrated  by  numerical  experiments  and  observation.  Numerical 
experiments  conducted  on  2D  models  with  embedded  scatterers  appear  to  explain 
coda  development  and  some  of  the  properties  of  the  Lg  phase. 

In  .section  2,  we  calculate  surface  wave  magnitudes  for  1 12  Nevada  Test 
Site  (NTS)  explosions  from  a  data  set  of  regional  long-period  seismograms  from 

North  American  stations.  In  order  to  utilize  the  nearer  regional  stations  (A<250),  a 
new  method  for  determining  Ms  and  been  developed  which  employs  synthetic 
seismograms  to  establish  a  relationship  between  the  amplitude  of  the  regional  Airy 
phase,  or  Rayleigh  pulse,  of  the  data  and  an  associated  surface  wave  magnitude, 
based  on  conventional  Ms  determinations,  calculated  form  a  synthetic  seismogram 
propagated  to  40°.  This  method  easily  lends  itself  to  implementing  path 
corrections.  The  inclusion  of  path  corrections  decreases  the  Ms  variance  by  a  factor 
of  two  and  also  decreases  the  average  value  by  0.08  magnitude  units.  This  latter 
effect  is  attributed  to  the  particular  station  network  used.  The  method  gives  stable 
Ms  values  that  correlate  well  with  other  magnitude  scale  values  over  a  range  of  three 
orders  of  magnitude  in  source  yield.  Our  most  refined  Ms  values  yield  the 
relationship  mb  =  0.9  x  Ms  +  1.82.  When  events  are  grouped  with  respect  to 
source  region,  significantly  better  fits  to  these  individual  site  linear  regression 
curves  are  obtained  compared  to  the  fits  obtained  using  a  single,  all  inclusive 
model.  This  observation  implies  that  shot  site  parameters  and  source  structure 
effect  surface  wave  magnitude  measurements,  although  event  yield  site  distribution 
al.so  may  be  in  part  responsible. 

Since  our  magnitude  values  are  based  on  a  theoretical  continental  structure, 
we  regres.sed  our  values  with  the  more  standard  values  of  Marshall  et  al.  (1979). 
Using  14  common  NTS  events  we  found  that  our  values  were  greater  by  0.53  +  - 
0.03  magnitude  units. 

In  section  3,  we  present  a  note  on  the  relocation  of  Tibetan  earthquakes. 
Present  ISC  legations  for  earthquakes  beneath  Tibet  indicate  a  random  distribution 
of  events  down  to  a  depth  of  al^ut  .50  km.  This  distribution  might  be  expected 
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from  a  relatively  cool  crust  which  would  allow  the  seismo-genic  zone  to  extend  to 
such  depths.  A  detailed  investigation  of  the  Tibetan  earthquakes,  with  magnitude 
greater  than  5.5  from  1964  to  1986,  yields  a  distinctly  different  picture.  Waveform 
modeling  of  depth  phases  such  as  ^  indicates  that  only  three  or  four  events  from 
this  population  is  actually  deeirer  man  25  km.  These  few  events  occur  near  the 
edges  of  the  Plateau  where  active  subduction  is  occurring  as  suggested  by  the 
thrust-like  nature  of  their  mechanisms.  The  events,  averaging  me  entire  population, 
occurred  earlier  than  indicated  by  me  ISC  by  about  3  seconds  which  leads  to  about 
a  1 .5%  and  0.5%  over  estimation  of  Pn  and  Sn  velocities  respectively  applying  ISC 
tables  and  standard  flat-layered  models.  A  more  serious  error  occurs  if  the  Pn  and 
Sn  velocities  are  determined  by  correcting  for  source  depm  but  assuming  me  ISC 
origin  times. 
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Broad-band  modeling  of  regional  seismograms; 
Imperial  Valley  to  Pasadena 

D.  Helmberger,  R.  Stead,  Phyllis  Ho-Liu  and  D.  Dreger 


Abstract 

This  paper  addresses  broad-band  wave  propagation  along  a 
corridor  from  Imperial  Valley  to  Pasadena,  California.  The  path 
consists  of  50  km  of  slow  basin  structure  with  a  shallow  moho 
followed  by  250  km  of  relatively  normal  Southern  California 
structure  approaching  Pasadena.  Events  occurring  in  the  valley 
produce  extended  body  wave  codas  which  carry  a  dispersive 
waveform  imprint  along  with  longer  period  Rayleigh  waves.  Many 
features  of  these  seismograms  at  periods  greater  than  a  few  seconds 
can  be  modeled  by  applying  a  finite-difference  technique  to  a  2D 
structure.  Shorter  periods  display  complex  behaviors  not  easily 
modeled  by  present  techniques,  especially  for  shallow  events. 

Events  occurring  in  the  normal  structure  produce  broad-band  SH 
waveforms  that  can  be  modeled,  analytically,  if  we  include  some 
extra  empirically  derived  parameters  determined  by  calibration 
events.  These  parameters  are  determined  by  adjusting  travel  time 
differentials  between  the  decomposed  wavefield  S,  sS,  etc.  For 
example,  delays  in  sS  relative  to  S  correspond  to  local  corrections  for 
sediment  cover  for  a  particular  source  region.  Asperities  and 
directivity  become  increasingly  important  at  shorter  periods,  as 
demonstrated  by  numerical  experiments  and  observation.  Numerical 
experiments  conducted  on  2D  models  with  embedded  scatterers 
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appear  to  explain  coda  development  and  some  of  the  properties  of 
the  Lg  phase. 


Introduction 

One  of  the  fundamental  reasons  why  quantitative  seismology 
and  waveform  modeling  is  contributing  to  source  retrieval  and  to  the 
monitoring  of  nuclear  weapons  is  its  prediction  capability.  That  is, 
given  a  set  of  seismograms  we  can  correct  for  propagational  effects 
and  predict  the  nature  of  the  source.  This  is  relatively  easy  when 
dealing  with  teleseismic  body  waves  and  well  dispersed  long  period 
surface  waves.  But  since  small  events  cannot  be  seen  teleseismically, 
we  must  address  the  more  difficult  problem  of  source  retrieval  from 
regional  phases. 

An  intermediate  stage  in  this  development  is  to  examine  broad 
band  records,  BB,  at  regional  distances  for  events  large  enough  to  be 
well-recorded  teleseismically.  This  feature  allows  the  source 
parameters  to  be  determined  by  modeling  teleseismic  bodywave 
phases,  direct  P  and  reflected  phases  (e.g.  pP,  etc.).  Starting  with  a 
known  source  allows  detailed  studies  of  regional  seismograms  in 
terms  of  path  effects  and  the  development  of  BB  Green's  functions. 
Hopefully,  these  Green's  functions  can  be  used  to  estimate  smaller 
events  in  the  general  vicinity  of  the  larger  event  or  master  event. 
Such  a  procedure  has  been  used  earlier  to  calibrate  surface  wave 
paths,  Romanowicz  (1982).  Unfortunately,  there  are  not  many  BB 
recordings  of  significant  events  presently  available.  Some  results 
from  detailed  studies  of  these  existing  records  are  encouraging  while 
others  are  not. 
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An  example  of  the  former  is  given  in  a  paper  by  Zhao  and 
Helmberger  (1991),  see  figure  1,  for  a  path  from  Quebec  to  Harvard. 
Their  strategy  for  modeling  these  records  was  to  break  the 
seismograms  into  segments  Pnl,  containing  Pn,  pPn,  sPn,  PmP. 
coupled  PL  waves;  Snl,  containing  Sn,  sSn,  SmS,  (etc.);  and  the 
fundamental  Rayleigh  waves.  Synthetics  are  generated  for  the 
various  segments  with  generalized  rays,  normal  modes  and 
reflectivity  methods  where  the  advantages  of  each  technique  are 
exploited.  Information  about  the  upper  crust  is  obtained  from  the 
fundamental  modes  and  crustal  thickness  and  velocity  gradients  in 
the  mantle  from  Pnl  and  Snl  waves. 

It  is  rather  remarkable  that  a  flat  layered  model  (see  Table  1) 
can  explain  so  many  features.  In  fact,  a  simple  layer  over  a  half 
space  does  a  good  job  for  the  first  100  seconds  of  long  period 
motions. 

These  results  suggest  that  in  stable  regions  where  Sn  propagates 
to  large  distances,  the  possibility  of  using  sSn  to  help  fix  source 
depths  could  prove  useful.  Note  that  Snl  and  sSn  are  considerably 
larger  than  Pn  and  pPn  and  can  probably  be  seen  for  smaller  events. 
The  usefulness  of  the  interference  produce  by  Pn  and  pPn  has 
already  been  demonstrated  by  Burdick,  Saikia  and  Smith  (1991)  for 
explosions. 

Explaining  the  jitter  in  the  BB  observations  displayed  in  figure  1, 
(assuming  a  flat-layered  model),  has  proven  difficult.  Thus,  it 
appears  appropriate  at  this  stage  to  introduce  scattering  into  the 
synthetics  as  proposed  by  Flatte  and  Wu  (1988),  Kennett  (1989)  and 
others.  Note  that  this  method  supposes  that  a  first-order  model 
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already  exists.  We  will  show  that  establishing  these  models  is  no 
easy  task,  especially  in  tectonic  regions,  as  discussed  for  a  BB  dataset 
for  a  corridor  from  Imperial  Valley  to  Pasadena.  Complexity 
abounds  along  this  profile  and  the  short  period  modeling  results  are 
not  so  encouraging. 

Observations 

The  observations  addressed  in  this  study  consist  of  a  collection 
of  low-gain  short  period  records  from  moderate  sized  events  along 
the  San  Jacinto  fault  system  which  have  known  source  mechanisms, 
and  BB  observations  of  recent  smaller  events.  Figure  2  gives  the 
locations  of  these  events  which  range  (roughly)  from  170  to  270  km. 
The  more  distance  events  are  beneath  the  Imperial  Valley.  Table  2 
list  some  of  the  event  information,  source  depth,  etc,  These  events 
are  all  strike-slip  in  nature  and  have  been  studied  in  detail  by 
several  authors:  for  example  see  Bent  and  Helmberger  (1991). 

One  of  the  reasons  for  working  with  the  Pasadena  data  is  the 
presence  of  a  recently  installed  Streckeisen  instrument.  Example 
seismograms  from  this  system  are  displayed  in  figures  3,  4,  and  5. 
The  bottom  three  rows  of  simulations  correspond  to  three  older 
types  of  instruments  operated  at  Pasadena  and  bracket  the  range  of 
frequencies  of  interest.  These  three  events  are  located  along  the 
extended  San  Jacinto  fault  system  and  can  be  used  to  compare  with 
the  older  well-studied  large  events  such  as  the  events  located  in 
figure  2. 

Since  the  azimuth  towards  Pasadena  is  near  the  (P-SV)  node 
assuming  a  strike-slip  orientation,  we  would  expect  considerable 
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variation  on  the  vertical  component  which  is  apparently  observed, 
(see  figures  4  and  5),  where  the  P-waves  appear  to  be  of  opposite 
polarity.  The  tangential  motions  are  similar  as  expected  along  this 
azimuth.  Recordings  from  the  aftershock  of  the  Superstition  Hills 
event,  see  (c)  in  figure  2,  display  distinctly  different  properties  from 
the  first  two  sets.  Note  that  the  BB  tangential  motions  are  longer 
period  and  the  ratio  of  wa.lp-to-wa.sp  period  torsion  responses 
increase  from  (3)  to  (5).  Stead  (1989)  compares  these  ratios  for  20 
events  in  the  region  and  finds  this  same  change  on  average.  Another 
difference  is  the  enhanced  strengths  of  surface  waves  for  basin 
events  and  their  relatively  slow  velocity  relative  to  to  the 
bodywaves. 

Note  that  as  the  Rayleigh  waves  become  stronger,  the  SV  waves 
become  weaker.  These  same  general  features  can  be  seen  in  the 
older  records  as  displayed  in  figure  6.  The  wa.sp  records  are  from 
the  large  events  and  the  wa.lp  is  from  a  smaller  aftershock  as 
indicated  in  Table  2.  Comparing  the  wa.sp  waveforms  for  a  large 
event  with  a  wa.lp  waveforms  of  a  aftershock  or  smaller  event 
generally  shows  good  agreement  as  displayed  in  figure  7.  This 
similarity  is  produced  largely  by  the  differences  in  faulting  area 
where  the  distributed  nature  of  the  main  event  tends  to  produce  an 
integrated  effect.  This  subject  becomes  interesting  because  of 
directivity  and  the  sensitivity  of  the  wa.sp  waveforms  to  fault 
asperities.  We  will  address  these  issues  in  this  paper,  along  with  the 
development  of  depth  phases  and  the  effects  of  random  scatterers  in 
the  crustal  waveguide. 
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Modeling  events  in  the  Imperial  Valley 

We  have  not  had  much  success  in  modeling  Imperial  Valley 
events  with  uniform  layers,  see  Ho-Liu  and  Helmberger  (1989). 
Finite-difference  modeling,  however,  appears  to  explain  some  of  the 
observed  long  period  features,  see  Ho-Liu  and  Helmberger  (1989). 
Their  model  is  displayed  in  figure  8  and  is  based  on  studying  a 
profile  of  long-period  (30-90)  observations  along  this  same  corridor. 
This  model  allows  events  occurring  in  the  basin  to  have  prolonged 
wavetrains  at  the  Pasadena  station,  PAS  while  events  north  of  the 
boundary  to  have  rather  sharp  pulses,  see  figure  2. 

Synthetics  for  basin  structures  can  be  generated  by  finite- 
difference  methods  following  the  approach  discussed  by  Helmberger 
and  Vidale  (1988).  The  basic  technique  is  based  on  the  expansion  of 
the  complete  three-dimensional  solution  of  the  wave  equation  in 
cylindrical  coordinates  in  an  asymptotic  form  which  provides  for  the 
separation  of  the  motions  into  SH  and  P-SV  systems.  Closed  form 
expressions  appropriate  for  finite-difference  source  excitation  are 
obtained  from  two-dimensional  Cagniard-de  Hoop  theory  with  rather 
elaborate  near-field  terms,  see  Vidale  and  Helmberger  (1987).  Note 
that  the  synthetics  generated  by  this  mapping  preserve  the  (P-SV) 
and  (SH)  systems  or  two  dimension  scattering  but  fail  to  treat  the 
SV-SH  scattering  problem  which  becomes  increasingly  important  at 
the  higher  frequencies,  see  figures  3,  4  and  5. 

Profiles  of  synthetics  for  the  geometry  displayed  in  figure  8  are 
given  in  figure  9  for  the  three  pass-bands.  The  first  fifty  km  of 
propagation  in  the  slow  Imperial  Valley  structure  causes  the  rapid 
development  of  dispersion.  After  the  shallow  Lovewave  leaves  the 
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basin,  it  speeds  up  and  disperses  more  slowly.  Three  component 
synthetics  generated  for  a  source  depth  of  12  km  appropriate  for  the 
1/28/88  event  are  given  in  figure  10.  The  overall  timing  looks  good 
but  the  shorter-periods  are  not  matched  well,  especially  on  the 
tangential  component. 

Comparing  these  synthetics  with  the  records  of  event  1  of  the 
1987  sequence,  trace  4  of  figure  6,  displays  about  the  same  degree  of 
fits  as  discussed  earlier,  see  Bent  ct  al.(1989)  where  this  event  was 
studied  in  detail. 

Results  for  the  Westmoreland  event  are  given  in  figure  11. 

These  results  are  not  bad,  but  the  short  period  complexity  becomes 
very  severe  for  shallower  events,  as  displayed  in  figure  12.  Some 
progress  in  explaining  the  long  period  motion  for  very  shallow 
events  is  given  in  Ho-Liu  and  Helmberger  (1989)  where  very  slow 
surface  waves  in  the  sediments  can  be  scattered  into  the  coda  by 
valley  edge  effects.  However,  generating  BB  synthetics  for  these 
types  of  models  is  very  expensive  relative  to  analytical  techniques 
and  therefore,  we  will  address  events  occurring  to  the  North  of  the 
Imperial  Valley  in  the  remaining  sections  assuming  more 
conventional  layered  models. 

Waveform  modeling  north  of  the  Imperial  Valley. 

In  this  section,  we  will  discuss  the  analytical  decomposition  of 
the  wavefield  in  terms  of  generalized  rays,  see  Helmberger  (1983) 
for  a  discussion  of  the  method  employed.  This  method  allows 
separation  of  the  wavefield  into  energy  that  starts  downward 
relative  to  upward  and  isolation  of  depth  phases  such  as  sS.  The 
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latter  phase  (sS)  has  potential  for  source  discrimination  since 
explosions  are  obviously  shallow. 

First,  we  examine  the  behavior  of  a  smooth  model  followed  by 
that  of  a  coarsely  layered  model  with  fewer  parameters  and  discuss 
their  properties  when  dipping  layers  are  allowed.  A  particularly 
effective  wavefield  decomposition  is  presented  in  terms  of  travel 
paths:  namely, 

1.  Direct  arrival  plus  surface  layer  multiples  (shallow  Love 
waves) 

2.  Diving  energy  paths  (lower  crustal  triplications) 

3.  Surface  reflected  paths  which  turn  below  the  source  (sS) 

These  latter  two  wave-packets  are  particularly  affected  by  vertical 
directivity,  as  demonstrated  by  broadband  experiments,  and  behave 
as  if  they  had  different  time  histories.  In  short,  combining  the  three 
contributions  with  various  time  shifts  which  correspond  to  depth 
changes  or  lateral  variation  allows  some  extra  freedom.  Modeling 
observations  with  these  added  parameters  becomes  much  easier  as 
will  be  demonstrated  later. 

A  velocity  model  composed  of  approximately  50  layers  and 
obtained  by  smoothing  a  southern  California  model  proposed  by 
Hadley  and  Kanamoii  (1979),  is  given  in  figure  13.  Figure  13  also 
displays  two  generalized  ray  sets  used  in  constructing  the  wavefield, 
namely  the  down-going  ray  set  and  up-going  ray  set  (excluding  the 
direct  arrival).  The  upper  portion  of  figure  14  displays  the  various 
contributions  of  these  three  raysets  to  the  total  potential  field.  These 
responses  were  produced  by  applying  the  Cagniard-de  Hoop 
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technique  to  the  generalized  rays,  see  Helmberger  and  Malone 
(1975).  The  upper  row  shows  the  loss  of  short-period  energy  with 
range  as  the  signal  becomes  diffracted.  The  down-going  rays  or 
diving  rays  contribute  significantly  to  the  high-frequency  content. 

The  moho  reflection  and  head  wave  complicate  the  picture,  especially 
due  to  contributions  from  sS.  The  interplay  between  S  and  sS 
changes  dramatically  when  we  change  mechanism  to  a  dip-slip 
because  of  the  polarity  shift  across  the  horizontal  axis. 

This  feature  is  easily  seen  in  some  preliminary  dipping  models 
used  in  modeling  waveform  data  from  the  1969  and  1968  events, 
see  figures  15  and  16.  In  this  situation,  we  decomposed  the  ray  set 
into  two  groups,  those  containing  the  direct  ray  plus  multiples  in  the 
surface  layer  (Lovewave)  and  those  diving  below  the  source  either 
directly  or  following  a  bounce  off  the  free  surface,  see  figure  15.  The 
direct  arrival  and  diving  rays  have  the  same  polarity  for  strike-slip 
events  but  opposite  in  the  dip-slip  case.  For  this  reason  the  dip-slip 
case  produces  much  different  looking  synthetics  as  displayed  in 
figure  15. 

A  comparison  of  synthetics  generated  from  this  dipping  model 
are  presented  in  figure  16  where  the  depth  effects  are  emphasized 
by  the  ray  decomposition.  Note  that  in  the  Coyote  Mtn.  example,  the 
source  depth  is  14  km  and  the  shallow  Lovewaves  become  less 
important.  For  the  shallower  Borrego  Mtn.  event  we  see  a  stronger 
Lovewave  contribution.  The  sSmS  phase  is  especially  obvious  in  the 
Coyote  Mtn.  record. 

One  of  the  advantages  in  working  with  large  events  in 
developing  a  model  is  the  constraints  on  source  parameters  obtained 
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from  the  teleseismic  modeling  of  body  waves.  The  Coyote  Mtn.  event 
appears  very  simple  teleseismically  and  occurred  at  the  relatively 
deep  depth  of  14  km,  see  Bent  and  Helmberger  (1991a).  A 
disadvantage  of  large  events  is  seen  in  the  Borrego  Mtn.  case  where 
the  source  contains  at  least  two  spikes.  This  complexity  is  clear  in 
the  teleseismic  results  as  well,  see  Ebel  and  Helmberger  (1982). 

The  dipping  model  discussed  above  was  derived  in  conjunction 
with  a  cross-section  running  from  Imperial  Valley  to  Pasadena  as 
displayed  in  figure  8.  The  model  predicts  synthetic  that  fit  these 
waveforms  quite  well,  although  it  is  poorly  constrained.  In  fact, 
there  are  only  three  principal  pulses  controlling  these  waveforms  at 
this  range.  They  are  the  shallow  Lovewave,  SmS  and  sSmS.  Since  a 
flat-layered  model  has  a  broader  applicability,  we  will  address  such 
a  model  next. 

Unfortunately,  we  do  not  have  a  profile  of  shear  wave 
observations  from  a  calibrated  event  and,  thus,  we  must  constrain 
the  model  from  other  sources  and  keep  the  parameters  to  a 
minimum.  Our  preferred  model  is  given  in  Table  3  and  has  just  five 
layers,  similar  to  Hadley  and  Kanamori  (1979).  Here  the  gradient  at 
the  top  of  the  model  given  in  figure  13  was  replaced  by  a  single 
layer  with  a  thickness  of  4  km.  This  thickness  was  adopted  from  a 
study  of  local  events  recorded  at  PAS  where  a  single-layered  model 
fits  the  three  component  data  very  well,  see  Dreger  and  Helmberger 
(1990).  The  smooth  positive  gradient  at  mid-depths  displayed  in 
figure  13  was  replaced  by  a  single  layer  with  velocity  3.6.  This  is 
consistent  with  studies  from  NTS  to  Tucson;  see  for  example, 

Langston  and  Helmberger  (1974).  The  transition  zone  at  the  base  of 
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the  crust  was  added  to  reduce  the  strength  of  the  moho  triplication 
as  suggested  by  the  P-wave  data,  (Hadley  and  Kanamori  (1979).) 

The  problem  discovered  in  modeling  the  waveforms  is  that  of 
preventing  the  moho  triplication  from  completely  overwhelming  the 
relative  weak  arrival  from  the  nearly  constant  velocity  crust. 

Applying  the  above  constraints  and  varying  the  thickness  to 
match  the  timing  and  waveshapes  produced  the  model  given  in  Table 
3.  Synthetics  appropriate  for  this  model  are  displayed  in  figure  17 
along  with  a  decomposition  which  proves  useful  in  modeling  near- 
regional  data.  The  first-column  shows  the  development  of  the  Love 
wave  associated  with  the  surface  layer.  Note  that  it  also  contains  the 
direct  arrival.  The  second  column  gives  the  results  contributed  by 
the  down-going  arrivals  and  the  third  displays  the  build  up  of  sS.  It 
is  clear  that  the  direct  wave  plus  the  surface  multiples  dominate  the 
motions  at  epicentral  distances  less  than  100  km.  This  feature  was 
observed  and  modeled  earlier  by  Helmberger  and  Malone  (1975)  for 
the  Morgan  Hill  earthquake  in  central  California. 

Most  of  the  San  Jacinto  events  occurred  near  distances  at  200 
km  where  all  three  raysets  contribute,  sometime  destructively  and 
sometime  constructively.  Adding  random  velocity  anomalies  to  the 
various  crustal  layers  produces  shifts  in  timing  between  these 
groups,  which  in  turn  yields  dramatic  high  frequency  effects  as 
discussed  in  the  next  section. 

Changing  the  source  depth  also  has  a  strong  effect  on  the 
excitation  of  Lovewaves  and  on  the  separation  in  arrival  times 
between  S  and  sS.  Figure  18  displays  the  depth  sensitivity 
appropriate  for  a  distance  of  205  km,  where  we  show  BB,  long-period 
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and  short  period  Wood-Anderson  synthetics.  The  phase  sS  is  the 
most  obvious  for  deep  events  where  the  Lovewaves  are  the  weakest 
and  less  interference  occurs. 

Application  of  the  decomposed  wavefield  with  "time  shifts"  is 
presented  in  figure  19.  The  fits  to  the  data  were  obtained  by  shifting 
the  phase  sS  back  in  time  by  .3  seconds  for  the  7/02/88  event  and  .5 
seconds  for  the  5/17/88  event.  Although  these  shifts  were  applied 
as  rather  arbitrary  corrections,  they  can  be  pictured  as  upper  crustal 
adjustments  for  these  particular  locations.  For  instance,  Hamilton 
(1970),  reports  about  .4  km  of  sediments  in  the  Borrego  Spring 
Valley  near  the  epicenter  of  the  5/17/88  event.  Thus,  sS  would  be 
delayed  about  .5  seconds  relative  to  the  diving  energy  in  accordance 
with  the  .5  sec.  lag  applied.  These  two  events  were  sufficiently  small 
to  treat  as  point  sources  and  apply  a  simple  trapezoid  time  history, 
namely  6ti,6t2  and  6t3  yielding  .2,  .4,  and  .2  seconds  for  both 
events.  Their  moment  estimates  are  4.5x10^1  ergs  for  the  7/02/88 
event  and  sightly  smaller  for  the  5/17/88  earthquake,  4.5x1021  ergs- 
Since  we  have,  presently,  only  one  broad-band  instrument  we  must 
treat  these  moment  results  as  preliminary  estimates. 

For  larger  events,  we  expect  considerably  more  complications 
involving  both  directivity  and  a  slip  distribution  containing 
asperities.  Figure  20  displays  predicted  synthetics  for  a  number  of 
rupture  geometries.  A  point  source  synthetic  is  included  for 
comparison.  The  case  of  forward  rupturing  event  is  nearly  the  same 
as  the  point  source  if  the  rupture  velocity  is  near  the  body  velocity 
as  usually  assumed.  Rupturing  away  from  the  receiver  produces  a 
strong  reduction  in  short  period  signals,  as  expected.  Downward 
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rupture  tends  to  enhance  the  moho  reflection  while  rupture  upward 
strengthens  sS.  Since  both  the  waveshapes  and  relative  short  period 
to  long  period  amplitude  ratios  are  so  sensitive  to  rupture  direction 
it  would  appear  that  broad-band  data  will  prove  particularly  useful 
for  the  purpose  of  studying  rupture  direction. 

If  we  scale  the  motions  given  in  figure  20  to  those  appropriate 
for  the  69  event  discussed  earlier  with  a  Mo=4.8xl0^‘^  from  Bent 
and  Helmberger  (1991a)  we  obtain  an  estimate  of  168  cm  for  the 
case  of  upward  rupture.  The  observed  peak  motions  were  53  cm  and 
the  waveform  fit  is  quite  good,  see  figure  21.  A  .30%  strength  of 
short  period  energy  release  to  long  period  level  is  about  the  average 
behavior  of  many  California  earthquakes,  (Bent  and  Helmberger 
(1991b)).  However,  to  avoid  biasing  this  estimate  by  directivity,  we 
clearly  require  multiple  samples  in  different  locations,  as 
demonstrated  in  figure  20.  This  will  be  possible  in  a  few  years  with 
modern  instrumentation,  at  least  in  California  with  the  installation  of 
Terrascope. 

Numerical  Scattering  Experiments 

The  introduction  of  laterally  varying  structure  causes 
considerable  interference  between  the  three  ray-groups  discussed 
earlier.  This  is  especially  true  at  the  larger  ranges  where  the 
summation  has  about  the  same  amplitude  as  each  individual  group. 

In  this  section,  we  discuss  the  theoretical  results  from  models 
containing  random  scatterers  embedded  in  the  various  layers.  We 
are  particularly  concerned  with  the  development  of  high  frequency 
coda,  and  with  amplitude  decay  with  distance. 
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The  model  and  geometrical  setup  is  displayed  in  figure  22.  We 
investigated  four  situations  with  variation  allowed  in  particular 
layers;  namely  RMOHO  (thin  layer  at  the  top  of  the  mantle),  RSURF 
(surface  layer  only),  RBOTH  (bottom  transition  layer)  and  RALL 
(variations  in  all  the  layers).  The  scattering  model  follows  the 
scheme  discussed  by  Frankel  and  Clayton  (1984)  where  a  Gaussian 
correlation  function  with  correlation  distances  (X)  in  the  vertical  and 
(Y)  in  the  horizontal  dimensions  is  applied.  Variations  in  velocity  are 
up  to  20%.  Locations  of  variation  are  random  and  their  positions  are 
indicated  in  the  above  figure.  Variations  of  this  size  have  been 
suggested  by  Frankel  and  Clayton  (1984)  and  others  and  are 
probably  on  the  high  side. 

The  complete  wave-field  was  generated  for  these  models 
ranging  in  distance  from  50  to  275  km  at  intervals  of  5  km  with 
frequencies  up  to  about  2  hz.  We  will  concentrate  primarily  on  the 
cross-over  ranges  as  displayed  in  figure  23  for  the  SH  field.  The 
synthetics  displayed  on  the  left  are  difficult  to  distinguish  from  the 
synthetics  discussed  earlier  from  the  homogeneous  model,  see  figure 
17.  The  next  three  columns  of  the  above  figure  shows  the  effects  of 
adding  random  scatterers  to  the  various  crustal  layers.  Scatterers  in 
the  top  layer  show  a  much  stronger  effect  than  do  those  in  the 
bottom  crustal  layer  (3.75  km/sec)  since  the  second  and  third 
columns  are  nearly  the  same.  Adding  in  scatterers  in  the  main 
crustal  layer  (3.6  km/sec)  causes  the  strongest  distortions.  The 
relative  amplitude  variation  across  the  various  rows  reflects 
primarily  the  velocity  at  the  receiver  which  is  the  same  for  the  three 
columns  on  the  right. 
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It  becomes  difficult  to  interpret  these  complex  waveforms  since 
we  can  no  longer  decompose  the  waveform  into  subgroups,  though 
we  can  still  identify  some  of  the  more  important  phases  such  as  SmS 
and  sSmS.  A  brief  review  of  figure  17  indicates  that  sSmS  becomes 
strong  near  200  km  and  SmS  and  sSmS  are  particularly  obvious  at 
the  ranges  170  and  185  in  the  first  three  columns.  Direct  SmS 
becomes  less  strong  in  the  most  severe  case  given  on  the  right. 
Apparently  as  the  ray  paths  flatten,  they  become  more  sensitive  to 
lateral  variation  as  one  would  expect  and  the  multiples  which  travel 
more  nearly  vertical  become  important. 

The  introduction  of  random  scatterers  into  the  surface  layer 
produced  significant  amplitude  anomalies  in  the  tangential  case  but 
did  not  greatly  increase  the  coda  which  is  so  commonly  observed  in 
regional  data.  The  effects  on  the  (P-SV)  system  are  more  dramatic  as 
displayed  in  figure  24.  Synthetics  for  the  RMOHO  case  are  displayed 
on  the  left  for  comparison.  Note  the  simplicity  of  the  P  arrivals, 
essentially  PmP  and  pPmP»  which  are  not  affected  much  by  the 
shallow  velocity  anomalies,  nor  are  the  corresponding  SV  phases. 

Significant  short  period  scattering  occurs  following  SV  which  is 
apparently  caused  by  the  interaction  of  the  Rayleigh  wave  with  the 
irregular  fine  structure  near  the  surface.  The  latter  contributes  to 
the  mysterious  Lg  phase  which  is  prevalent  on  short  period  regional 
records.  Convolving  the  broad-band  SH  responses  of  figure  23  with  a 
(wa.sp)  instrument  yields  similar  looking  records,  but  showing 
considerable  modulation  in  amplitudes,  as  displayed  in  figure  25. 

The  amplitude  fall-off  or  attenuation  as  determined  by  the 
strong-motion  community  reflect'  the  choice  of  crustal  model  The 
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thick  mid-cnistal  laver  of  velocity  (3.6  km/sec)  i.s  probablv  qu 'e 
common  in  western  United  States  and  is  responsible  for  the  rapid 
decay  between  50  to  90  km.  The  increase  at  larger  ranges  is 
influenced  by  the  lower  crustal  transition  and  the  sharpness  of  the 
moho.  Note  that  these  amplitudes  are  also  influenced  by  radiation 
pattern.  These  plots  are  appropriate  for  the  strike-slip  case  which  is 
the  most  common  type  observed  along  the  San  Jacinto  fault  system. 

The  strong  arrival  at  the  range  of  105  km  in  the  full-scattering 
model  is  caused  by  constructive  interference  between  the  moho 
reflection  and  first  multiple  in  the  top  layer  (Lovewave).  In  general, 
the  strong  scattering  introduced  in  this  exercise  was  sufficient  to 
obscure  the  moho  reflection,  but  the  possibility  of  large  motions  near 
the  moho-cross  over  is  hard  to  avoid. 

The  increased  amplitudes  near  100  km,  as  discussed  here,  are 
not  apparent  in  most  strong-motion  datasets.  However,  this  may 
reflect  the  prevalent  processing  methods  normally  applied;  where 
one  truncates  the  range  of  interest  at  the  first  strong-motion  station 
that  does  not  trigger.  In  light  of  the  strong  motions  observed  in  San 
Francisco  for  the  Loma  Prieta  (1989)  earthquake,  (Somerville  and 
Yoshimura,  1990)  and  observations  from  the  Saguenay  earthquake 
(1988),  Somerville  et.al  (1990),  this  subject  should  probably  be 
revisited. 

Conclusions 

This  paper  investigates  BB  wave  propagation  along  a  corridor  in 
Southern  California  from  the  Imperial  Valley  to  Pasadena. 
Seismograms  from  Imperial  Valley  events  are  characterized  by  SH 
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waves  which  appear  dispersed,  and  possess  well-developed  codas. 
Rayleigh  waves  are  dispersed  and  relatively  strong  compared  to 
body  phases.  These  features  can  be  modeled  with  a  finite-difference 
code  at  periods  greater  than  a  few  seconds  by  assuming  the  following 
two-dimensional  model;  a  slow  basin  with  a  shallow  moho  followed 
by  a  thick  crust.  Propagating  50  km  in  the  basin  is  sufficient  to 
disperse  the  wave-field.  Shallow  sources  apparently  excite  the 
surface  waveguide  where  the  field  travels  slowly  to  the  edge  of  the 
valley.  There  it  is  re-radiated  very  strong  coda  signatures.  It 
appears  difficult  to  propagate  the  wave-field  across  boundaries  of 
this  type  by  present  analytical  techniques. 

Modeling  events  at  smaller  ranges  in  the  normal  crust  proved 
possible  with  analytical  codes.  In  these  codes,  a  decomposition  of  the 
wavefield  into  the  following  individual  phases  was  useful:  direct  S 
and  Love,  diving  S  and  sS  in  the  tangential  case.  Shifts  in  SmS 
relative  to  sSmS  of  a  few  tenths  of  a  second  allow  excellent  fits  in  BB 
records  in  many  situations.  Directivity  plays  an  important  role  in  the 
sharpness  of  these  phases  which  adds  to  the  richness  in  modeling  BB 
regional  phases.  At  these  ranges  the  Moho  critical  angle  reflections 
contribute  the  strongest  peak  amplitudes,  at  least  in  flat-layered 
models.  Numerical  experiments  conducted  with  the  inclusion  of 
scatterers  into  the  various  layers  did  not  change  this  result. 

Scatterers  in  the  surface  layer  produce  signals  which  resemble 
typical  Lg  phases,  at  least  as  observed  on  the  vertical  component 
waveform. 

In  conclusion,  it  appears  possible  to  model  depth  phases  such  as 
sS  in  tectonic  regions  with  analytical  methods  if  major  geological 
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boundaries  are  not  crossed.  Secondly,  the  ratio  of  S  to  Love  waves 
can  be  used  to  estimate  source  depths  if  paths  are  calibrated.  Future 
efforts  will  include  the  application  of  source  inversion  techniques  to 
these  type  of  records. 
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Table  1.  Crustal  model  MPM 


Table  2.  Earthquakes  along  profile 


date 

origin  dme 
GCT 

latitude 
degrees  N 

longitude 
degrees  E 

depth 

km 

magnitude  distance 
from  Pas. 

03/23/54 

04.14 

33.28 

-116.18 

16.0 

5.1 

230  km 

04/09/68 

02.28 

33.18 

•116.12 

11.0 

6.4 

224 

04/28/69 

23.20 

33.33 

-116.33 

20.0 

6.1 

191 

04/26/81 

12.09 

33.13 

-115.65 

6.0 

5.3 

260 

1 1/24/87 

01.54 

33.08 

•115.78 

5.0 

5.7 

251 

11/24/87 

13.15 

33.01 

-115.84 

2.0 

6.0 

252 

01/28/88 

02.54 

32.91 

-115.68 

6.0 

4.6 

270 

05/17/88 

19.38 

33.24 

-116.25 

8.0 

3.8 

205 

07/02/88 

00.26 

33.49 

-116.44 

12.0 

4.1 

176 

03/06/89 

22.16 

33.17 

-115.59 

1.0 

4.7 

262 

Table  3.  Normal  Southern  Cal  Crust 


Figure  1.  The  top  pair  of  traces  compares  a  synthetic  with  the  first  100  sec  of  the  venical 
observation  of  the  Quebec  earthquake,  1 1/25/88,  as  recorded  at  the  Harvard  broad-band  station,  A 
=  625  km.  The  lower  three  traces  display  the  fits  after  convolving  with  the  WWSSN  long  period 
response.  The  first  30  to  40  secs  is  composed  mostly  of  P  diffracted  energy  along  the  top  of  the 
mantle,  Pnl.  The  latter  30  secs  is  composed  mostly  of  SV  diffracted  energy  along  the  same  mantle 
path,  Snl- 


Figure  2.  This  map  displays  the  locations  of  a  number  of  events  discussed  in  this  paper  relative  to 
the  extended  San  Jacinto  fault  system.  The  stars  indicate  relatively  large  events  while  the  circles 
indicates  smaller  recent  events  recorded  by  the  new  IRIS  system  at  Pasadena. 


Event:  07/02/88  00:26 
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'  30  s 

Figure  3.  Display  of  motions  recorded  from  the  7/2/88  magnitude  4  event  (a),  A  =  176  km.  Peak 
amplitudes  are  given  in  cm.  The  gain  of  the  Press-Ewing  has  been  set  at  one  in  these  simulations. 


Event:  05/17/88  19:38 
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Figure  4.  Display  similar  to  that  in  figure  4  of  the  5/1 7/88  event  (b),  A  =  205  km. 


Event:  03/06/89  22.16 
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'  30  s  ' 

Figure  5.  Di.splay  similar  to  that  in  figure  4  of  the  Supierstition  Hills  aftershock  1/28/88  event  (c), 
A  =  256  km. 


3/23/54  wa.lp 
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Figure  7.  Comparison  of  tangential  motions  of  large  events  (wa.sp)  with  small  events  (wa.lp). 


31 


Figure  9.  Synthetics  appropriate  for  a  strike-slip  source  beneath  Imperial  valley  (d  =  7  km) 
assuming  the  model  defined  in  figure  8. 
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f  10.  Comparisci.  r'T  synthetics  and  obsci^'ations  for  the  Superstition  Kills  aftershoc'k 
constructed  from  the  model  displayed  in  figure  8  assuming  a=1.73p. 


Figure  11.  WA  (lOOx)  recording  from  the  Westmorland  eanhquake.  Two  sources  were  used  in 
generating  the  synthetics,  one  with  a  M©  of  2.5  x  10^^  dyne-cm  at  a  depth  of  7  km  and  the  second 
one  with  an  Mo  of  4.5  x  I02^  dyne-cm  at  a  depth  of  10.5  km,  delayed  by  about  2.5  secs.  This 
2.5  sec.  delay  is  also  observed  in  local  strong  motion  accelerograms.  The  lower  comparisons  are 
made  after  filtering  the  synthetics  and  data  by  one  and  two  second  triangles. 
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Evenl;  03/06/89  22:16 

Tangential  Radial  Vertical 

1  ...  0.0021  1  n  nnri'y 
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F/gTO  12.  Display  of  motions  reconled  from  the  3/OT9  event,  (d  in  figure  2)  which  occurred  near 
the  Westmorland  event. 


Velocity  (km/s) 


f  igure  1 4.  Wave-field  decomposition  displaying  the  response  of  the  direct  arrival  on  the  top 
followed  by  the  contribution  from  down-going  paths  (S)  and  up-going  paths  (sS).  The  bottom 
row  displays  the  synthetic  assuming  a  (.2,.2,.2)  trapezoidal  source. 
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ing  siep-responscs  for  a  dipping  layered  model 


Fiirure  16.  Synthetics  constructed  from  the  step-responses  (strike-sHp)  displayed  in  figure 
15  for  a  range  of  19 1  km  along  with  a  similar  construction  for  A=224  km.  Comparison.s 
with  the  data  is  given  in  the  bottom  row.  The  time  function  for  the  Coyote  Min.  event 
consisted  of  a  trapezoid  (.2,2, .2).  The  time  function  used  in  the  Borrego  Mtn.  event 
'  oi’sisteii  of  two  trapezoids  with  the  same  timing  where  the  second  is  three  times  the  first 
aiul  laged  .7  sees. 
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Figure  18.  Depth  sensitivity  sections  {A=224  km)  showing  the  domination  of  the  Love 
wave  at  shallow  source  depths.  The  solid  line  indicates  the  Sn  arrival  and  the  dashed  the 
moho  reflection  (SmS)  and  corresponding  surface  reflected  phase  (sSmS). 


Figure  19.  Comparison  of  observations  (solid)  and  synthetics  (dashed)  for  event  (a)  and 
(b),  see  figure  2,  the  lower  panels  display  the  wave-field  decomposition  assuming  a 
(.2,.4,.2)  trapezoidal  source. 
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Figure  20.  Directivity  effects  are  simulated  by  assuming  propagation  along  a  line  at  six 
events,  one  km  apart,  rupturing  at  3  km/sec.  The  events  are  all  strike-slip  with  moments  of 
l.xlO-3  ergs  and  have  durations  of  (.2,.2,.2). 


Figure  21.  Comparison  of  the  upward  propagation  case  with  the  observed  Coyote  Mtn 
event.  The  peak  amplitude  of  the  data  is  .19mm  corrected  for  instrument  gain  (low  gain 
l(K)x),  predicted  to  have  a  peak  amplitude  of  53  cm  if  recorded  on  the  standard  Wood- 
Anderson  instrument  (28(X)x). 
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Figure  22.  Diagram  displaying  a  crustal  model  with  various  layers  containing  velocity 
anomalies  (20%  extremes). 


FD  with  Random  Media  Reducing  V=3.90 


F  igure  23.  Profiles  of  BB  tangential  synthetics  through  the  cross-over  distances  of  50  to 
2(i)  kni.  The  column  on  the  left  (RMOHO)  is  nearly  identical  to  the  homogenences  case, 
see  figure  17.  The  next  three  columns  include  increasingly  more  layered  anomalies. 


46 


hi 


Figure  24.  Comparison  of  BB  and  shon  period  synthetics  for  essentially  the  homogeneous 
layered  model  with  the  model  containing  anomalies  in  the  top  layer.  Note  the  increased 
short  period  duration. 


Amplitudes  of  Random  Media  Finite  Difference  Runs 
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caused  by  critical  an^le  Muho  reflections. 


SECTION  2 

Deteniiining  M.s  Magnitudes  from  Regiiwial  N'l'S  Data 
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Determining  Ms  Magnitudes 
from  Regional  NTS  Data 

Brad  Woods,  David  Harkritlrr 

August  2,  1991 


Abstract 

VVo  hav('  calciilalod  surface  wave  magnitudes  for  102  Nevada  Test 
Site  (NTS)  underground  nuclear  explosions  from  a  data  set  of  regional 
long-p<'riod  seismograms  from  North  yXmericaii  stations.  In  order  to 
utilize  the  tu'arer  r<'gioual  stalioii.s  (A  <  2')°),  a  new  met  hod  for  det(>r- 
miuing  has  been  develop«>d  which  <*mploys  syiitlietii  .-ei.smograiiis 
toestaiilish  :i  relationship  l)etwe«>n  the  amplitiid('  of  tiu'  regional  Airy 
phas<',  or  Uayl<'igli  pulse,  of  the  data  and  an  associated  surfac  e  wave 
magnitude,  hasc'd  on  conventional  Ms  detc'rminat  iotis,  calculated  from 
a  synthetic  seismogram  [rropagatc'd  to  dO®.  'I'liis  im'tliod  easily  lends 
itself  to  im|)lementing  (rath  corrections.  'I'he  inclusion  of  (rat  h  correc- 
lions  dec  Teases  the  Ms  variance  Iry  a  factor  of  two  and  also  dec  icasc's 
the  avc'ragc'  value  hy  O.OH  magnitude  units.  'I'his  latter  cdlec  t  is  at 
trihuted  to  the  (larticular  station  ueUwork  usc'cl.  The-  method  gives 
stable  Ms  Viduc's  that  correlate  wcdl  with  otlic-r  magnitude  scale  wdiies 
over  a  range’  of  thrc’e  orcic’rs  of  magnitude’  iti  source’  yic’hl.  Our  most 
rc'finc’cl  Ms  value’s  yic'ld  the'  rc’latie)nshi|>  iiii,  -  tl.so  •  M.  •  11.  wli'  ie  II 
is  cle[)c’nclc’nt  u|)c)n  source  rc’gion  ami  shot  im’dium.  I  his  relationship 
holds  for  e’vents  of  all  si/c’s.  WIm’ii  c’vc’iits  arc’  gron|)C’il  with  ie  s[)e(  t 
to  source  rc’gion.  signilicaiitly  bc’ttc’r  (its  ter  t  lu’sc’  indivieiual  siic  lin- 
c’ar  rc’grc’ssion  curves  are’ eebtaiiied  ce>m|»arc’d  to  the  (it s  editained  using 
a  single’,  all  inclusive’  menle’l.  'This  ob.sc’rvation  im|ilies  that  shot  site’ 
(raratnc’te’rs  and  seuirce’  structure’  elb’ct  siirlaee’  wave-  magnileiele  nieui 
sure’tne’nts,  although  eve’iil  yield  site’  distribution  also  may  be'  in  p.irt 
rc’,s|)onsil)|e’. 


■SO 


rlltirl  wciil  Uiwiirtls  iiuikiii)*  tln' sft  i'iiin|>r<'li)'iisiv<‘ 
<livf'is('  ill  Icrms  of  yii'ld,  soiirrc  location  and  sliol  iiifdiiiiii  in  order 
to  di'tenniiK'  tlie  portaliility  of  seismic  measuriiiK  scales.  In  particular 
we  exaniiiie  I'aliiile  .\Iesa,  Rainier  Me.sa  and  ^'ncca  F'lat  explosions 
deioiiated  aliove  and  lielow  llie  water  table. 

Since  oiir  ma^iiit iid<‘  values  are  based  on  a  tlu'oreliral  continental 
si  nil  lure,  we  regressi'd  our  values  with  the  more  standard  value's  of 
.\Iarsiiall  e/  a!  (Ifl7!))-  I  siiiff  Id  rommon  NTS  evenis  w<'  found  that 
our  values  w«'re  Rteaie'r  by  by  O.od-i--  0.03  magnitude  unils. 


1  Introduction 

\\('  re'-cxanilMc  I  lie  use  of  surfa<<‘-waves  for  ntl<^<'r!^rollIui  nuclear  I'xplosioii 
iiiagnit iide  deh'i  iiiinal  ions,  particularly  for  smaller  v  ii'ld  (N’  <  2()Kl)  evenis. 
I  lie  surface  wa\e  mattnil  ud<'-\ ie'ld  scaliiifi^  law  lor  sin  li  low  yield  e\enls.  until 
now,  was  not  known  well.  Tin'  <lata  ii.scel  an*  loiift  period  .Norlli  Amerii.in 
station  vertical  records  lor  102  spe'cifie'd  Nevada  I'esl  Site'  (N  I  S)  e'Ve'iits.  Tin' 
stations  u.se'el  arc  from  several  networks,  Tbeir  respective  iust  runteuls  all  l\ave' 
jrass  banels  that  lie'  within  the  10  to  60  second  range.  Snrfaee'  wave's  are'  ve'iy 
useful  for  yie'lel  e'stimation  purposes,  for  (M«)  is  de'te'rmine'el  from  re'lalive'ly 
long-period  .se'ismic  wave’s  whicli  are  insensitive  te)  high  fre'epie’iiry  ne'ar-soun  e' 
effects,  whie  h,  along  with  seve'ral  other  possible'  mechanisms,  may  lx*  cause'ei 
by  asynmu’t rie's  in  tlie  shot  cavity,  see  Zhao  and  Ilarkride'r  (1901).  I'he'se' 
high  fr<’c|uency  .serurce  effc'cts  may  cause  appreciable  bias  in  magnitude's  that 
are’  baseel  eui  higher  fre'cpie'iicy  waves,  such  as  the  im,  and  1.^  scales. 

For  the-  lower  yielel  e'vents  it  beconu's  ne’ce'ssary  te)  incluele'  the'  data  frerm 
ri'gional  statieins  (A  <  25°),  for  tele.se'ismic  re'eerreiings  have'  teio  leiw  a  signal 
to  noise  ratie).  which  make's  them  unusable.  At  re'gional  elistance's  snrfaee' 
wave’s  are-  not  we’ll  elisjie-rse’el,  having  a  prerminani  .Airy  phase’  pulse'  with  a 
pe'iieiel  be-lwe'e’n  6  aiiel  20  se’eernels  (Ale’wine,  1972).  se>  that  it  is  not  perssible' 
to  me'asnre’  .M,  e onve'iil ionally  (  that  is  me’asiiring  the’  am|>lil nele- of  the-  20 
see  .  wave').  Feir  Nejrth  .Ame'iiea  in  ge'u'e'ral,  tlie'ie'  is  mininum  in  the'  genip 
vi'locity  e  nrve’  ne-ar  12  .se’eernels  for  the  funelame’iilal  Hayle’igh-wave’  (  Marshall 
ft  III  1979). 

lei  measure’  M,  we  e-mplery  a  technie|ue  whereby  I  lu'ene’l  ie  al  se'ismograiiis 
in  e  einjunct  ion  with  I  he  surface  w’ave  data  arc  empleiye'el  to  iinlire’ctly  eahulale- 
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it.  Ill  u.sing  this  procedure  several  propagation  path  models  were  tested  to 
determine  the  effect  of  attenuation  and  seismic  velocity  structure  upon  the 
Ms  values.  Thesi'  calculated  M*  values  remain  stable,  have  r('asoiiat)l y  small 
errors  and  correlate  well  with  associated  mb  maginiliides  and  log  yield  for 
the  event  data  set.  The  M*  —  mb  relationships  are  determined  by  a  weighted 
least-.squariis  linear  regression. 

Besides  comparing  the  M*  results  with  several  different  independent  mag¬ 
nitude  scales,  the  data  have  also  been  separated  with  respect  toscaircc'  rc'gion 
and  shot  matcuial.  vahms  at  Yucca  Flat  lend  to  be  larger  than  those  at 
l{ain<'r  Mc'sa  by  ().()!•  magnitude*  units  for  a  given  mj,.  I  lu  re  also  a|j|>ears  to 
lie  seinu'  dillc'rc'iicc*  in  wavc'forins  between  evc'iil.s  of  these  two  source  regions. 
Fahiite  Mc’sa  events  aic*  0.22  magnitude*  units  large'i  than  those  at  'I’ucca  flat 
for  explosions  sc't-off  below'  the  water  table  and  with  the  same  nii,. 

We  do  not  account  for  t<*ctonic  reh'ase  effc'cts  upon  the  magnitude  mc'a- 
siirments.  Such  effects  are  bc'st  accounted  for  with  inornent  tensor  inversions 
of  sources  which  invoiws  more  .sophisticated  data  analysis.  .Standard  .\1, 
mc'asurement  techniciuc's  ignore  this  factor  as  wi'll. 

2  Data 

The  data  are  long-p<*rio(l  vc'rtical  seismograms  recorded  at  .North  Ami'iica 
stations  for  102  c'xplosions  at  NTS  and  consist  of  digitized  World  \\  ide  Sc'imic 
Network  (WWSN)  and  Canadian  Seisinographic  Network  (CSN)  records. 
Digital  World  Wide  Seismic  N(*twork  (l)WWSN),  bawrence  Livermore  He 
gioiial  Seismic  Network  (LNN)  and  |{c*gional  lest  Seismic  .Network  (HSI  .N ) 
digital  data  for  events  oi  c  urling  lat<*r  than  lOSI.  Th<*  analog  \\  W'S.N  and 
('anadian  station  data  w<*re  digitizc'd  by  KNSCO.  I'ifty  e  ight  stations  com¬ 
prise  the*  nc*tvvork,  although  fc'wc'i  than  00  percent  of  tin*  stations  had  data 
available  for  any  single  c'vc'iit.  Lig.  1  shows  a  map  of  the*  station  ncriwork. 
Fpicc’iit ral  distancc*s  range*  from  220  km  for  .N'I'S  to  (ISC  ( ( loldstone.  Cali 
fornia).  to  .'}20l)  km  for  N'I'S  to  S'lM  (Saint  .lohns,  .Newfoundlaiul ).  for  the 
smallc'r  c*vc’nts,  |)articularly  Hainer  Mc*sa  <*xplosions.  only  the*  nearc'r  stations 
(distance*  <  1000  km)  had  c*ithc*rdata  avaiiablc*  or  re  asonable  signal  to  noise- 
ratios.  Station  coverage  varies  widely  bc*twien  c-ve-nts.  l  ive*  of  the*  smaller 
c*vents  only  hael  one  vialrle  station  seismogram  e-ach,  while*  some*  e'vents  hael 
over  110.  I'lic*  a vc’ragc*  numbc*r  of  stations  reporting  pc'r  eve-nt  is  10. 
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S(irf«<  <'  wavc^  t  lial  ))r(jpaga1«‘  a<  r<)ss  orcanK  /coiil  iiiciilal  margins  ntulrrgo 
significant  modifical  ion  in  llicir  waveforms  Ix'cansc  of  I  lie  great,  lateral  vari¬ 
ation  in  crustal  aiul  iipix'r  maniK’  slnictiire  at  siicli  Iroimdaries.  'I'liese 
|)ro(>agat i' >11  elfects  are  not  straight  forwani  to  nuxh'l,  hence  approiniale 
(Ireen's  fnix  lions,  or  transfer  fuix  tioiis.  are  dillicnil  to  obtain.  Without  ro- 
liiist  (Iri-en's  functions  it  is  hard  to  infer  accurate  .source  information  from 
lh>'  data.  .Smaller  events  als<»  are  not  likely  to  he  ohservi-d  at  the  distant 
statiiiiis.  which  often  include  oceanic  st riictnre  along  tlx'ir  propagation  [tat  h, 
and  make  llx-sc'  longer  paths  c'vc'ii  le.ss  attractive*  to  include’  in  the  monitoring 
tie'twork.  Hence',  we*  chose'  to  confine*  onr  stiiely  to  surface'  wave’s  trave’lling 
soh'ly  alonii  eemt iiieiital  patlis.  i.e.  witliin  Nextli  .\me’rie'a. 

Of  the'  ll)i  e've'iits.  J?  are*  frexn  Pahiite  Mesa.  Hi  are*  freim  Haine’r  Me*sa,  .')S 
from  't’ue  e  a  f  lat  aitel  one*.  Pile'driver.  is  le)cate*el  at  (’limax  .Stoc  k.  We*  consiele*r 
the'se*  to  he*  I  distinct  sexiree*  re*gions.  I'br  some*  spc’cifie'  stations,  waveforms 
varied  soine’what  lx'twe*e*n  events,  depending  upon  soiire  e*  Icx  ation. 

I  ho  I’ilc'drivc'r  data  from  a  given  station  lex>k  airprc’ciahly  eliffere’nt  from 
that  of  any  othe'r  e’ve*nls  re’corele*el  at  that  same*  station.  I  his  was  true*  for 
c've'iv  station  ri'e  erreling  l’ile'elrive*r  and  |)rohahly  is  eanse’d  hy  iliire*re*nce*.s  in  the* 
source*  re’gion  for  this  e’xplosion.  Pile*elrive*r  was  de’toiiate’d  in  a  granitic  source* 
re'gion,  north  of  the*  othe’r  site's.  The’  source*  to  re*ee'ive*r  gc'ome’trie’s  for  this 
e’ve’nl  are’  a|)proxiiiial e'ly  I  lie’  same*  as  those*  as  the*  othe’r  N  I’.S  e’Ve'iits.  so  the’ 
difle’re'iice'  in  w  a  ve'forins  dex'sn  I  appear  to  he*  at  I  ri  lull  a  hie*  to  elis[>e’rsi\e*e*lfe’cls 
caiisc’el  h\  elilfe’re’iice’s  in  pro|>agation  path  length.  I’ilc'dri ve’r  was  the*  only 
Cliamx  Siex  k  e*vent  with  re'adily  available*  data,  .so  no  flirt he*r  e’xaniination  of 
tills  source’  was  e'arrie*el  out. 

.\t  se)ine’  of  the*  ne’are’i  re’gional  stations  (distanee*  <  !t°),  the*re*  are  also 
suhtle*  ditre’ie'iice’s  Ix’twe’e  n  the’  ^  uee a  I’lat  and  l{aine’r  Mc’sa  eve*iil  wave*forms. 
,\t  1)1(1  (Diigwav.  I'tah).  for  e’xample,  tin*  l?aine*r  e'vent  waveforms  lex>k  as 
if  the’  .\irv  phase  has  Ixe’ii  llillx’it  transforme*el  (e'ejuivale’iit  to  a  ht)®  phase* 
shift)  re’lative  to  the  Yucca  Mat  waveforms.  I'he  DUC  Kainer  wave’forms 
also  cejiitain  more*  high  fre*e|ue*ncy  cexia  energy  than  t  hose*  from  \’ue  e  a  I’  lat  . 
I’ahiite  e’ve’uts  are  similar  in  wave*form  to  Rainier  e*ve*nts  and  have*  le*ss  high 
fre’epie’ne  y  eexile’iil  than  Yucca  evc’iits. 

I  he*  se  ismograms  wi’ie*  hanel-[)a.sse*el  fillere*el  lx*twe*e*n  (i  and  10(1  se’coiiels 
III  gi’t  rill  III  the’  long  pe’lioil  aiiel  short  pe*riexl  iioise*  whiell  Wolllel  otlu’iwise 
aHe’ct  the*  pi'.ik  to  pe’ak  mi’asure’me’iit  of  the  Rayle*igh  pulse*.  I'he’  ah.soliile’ 
ami^lit iieh’s  of  the*  data  we’re*  verified  hv  chexrsing  se-veral  elifferent  stations 
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and  coinparing  tlu*  peak  to  peak  amplitude  of  the  digitiz<‘d  record  to  that 
mea,siire<l  <lire(  tly  from  tlie  resp«'ctive  analog  VV'WSSN  Him  (  hip.  I'rom  thi.s 
.step  we  ascertained  that  thestati(»n  gain  fa<  tor  had  heeri  (onci  lly  fac  tored 
out  of  the  sc'ismograms. 

3  Ms  Calculation  Technique 

We  have  developed  a  method  to  measure  surface  magnitiuh's  indirc'ct  1\'.  He 
(  au.se  a  large  |)ort  ion  of  the  data  for  low  yic'ld  ('vent  s  is  from  st  at  ions  recording 
at,  rt^gional  distanc«'s  (A  <  ‘25°)  ,  it  is  not  possible'  to  calculate  .\I,  conveii- 
tionally,  for  the  Hayh'igh  wav«'  is  pidse-lik**  which  prc'c  hides  measuring  a  well 
dispers<‘d  ‘20  se< .  pha.s<'  (Alwine,  197‘2).  WeaddK'ss  this  problem  with  the  use 
of  synthetic  seismograms  of  the  fundamental  Rayleigh- wave  using  a  modified 
llaskel-'Pliompson  matrix  method  (llarkrid<>r,  liKil). 

For  each  source  to  receiver  path  a  theoretical  Rayleigh  wave  is  generated. 
'I'he  Earth  modc'l  used  to  create  this  synthetic  is  meant  to  relh'c  l  the'  average 
I’atrth  st  ructure  between  N  I  S  and  the  given  station.  I'he  l-.aitli  models  isc'd 
in  this  study  wc’H’ determiiu'd  from  inversions  of  dispersion  and  allenualion 
data  as  well  as  forward  modeling  of  the  waveform  to  line  tune  the  models. 
I'he  criteria  for  determining  the  goodness  of  lit  of  the  synthelit  to  the  data 
are  dis|)ersion,  absolute  travel  time  and  waveform  fit  (relative  amplitude  of 
dilferent  dispersed  phases).  Hence  the  synthetic  seismogram  displays  the 
same  spectral  and  time  domain  waveform  characteristics  as  the  data  which 
it  simulates.  'I'his  was  done  f«u'  all  pat  hs.  The  pat  hs  to  \V\\  S.\  aiul  (^anadian 
stations  were  taken  from  a  study  by  Stevens  (Stevens.  IbSti).  1  lu'  RSI.N. 
bLN  and  DWWS.N  paths  were  determined  previously  by  this  research  grouj). 

'lb  determine  M^  for  a  particular  source- receivc'i  geometry  two  svnthetii  s 
are  g<'nerated.  One  which  is  propagated  the  a(tual  path  distaiKc  that  is 
meant  to  simulate  the  data  and  one  whi(  h  is  propagated  to  It)  .  .\t  10"  tin 
surface  wa\’e  train  is  well  dispersed  and  stable,  so  that  a  ( ()n\eni  i(»nal  .\b 
value  can  be  calculated.  To  calculate  M,.  w('  u.s<'  the  moditied  \  on  Seguern 
formula  (Vbn  Seggern.  1977): 

Ms  log,„(  \/l  I  i  i  I  d.^. 

where  A  is  the  peak  to  peak  amplitude  (in  nanometers)  of  the  wavelet  mea 
sur('(l  from  the  vertical  re(()rd,  1'  is  the  period  of  the  wavelet,  measured  in 


seconds,  and  A  is  the  {)ro|)agation  distance  in  degrees.  'I'his  .scale  was  chosen 
Ix'cause  I  lie  distance  coeliicie'nt  ( 1 .08)  more  clo.sely  approximates  tin*  (‘ffecl  of 
attenuation  along  contitK'tital  paths.  A  vertical  component  measurement  has 
two  advantages  over  liori/ontal  (<>m|>onent  iin'asnrements.  The  horizontal 
components  nsnallv  have  lovver  signal  to  noise  ratios  than  the  vertical  com 
poniTit  and  gi'iierally  aie  more  liki’lv  to  he  I'ontaiiiinated  hy  hove  vvavesignals 
which  mav  lie  generated  hy  ti'clonic  rel«*ase,  source  elfecls.  or  scattering  <lne 
to  lateral  variations  in  tlu'  Earth’s  structure. 

Both  till'  regional  and  teleseismic  synthetics  are  generated  with  the  same 
.soiiKi-  fiUK  tion.  so  that  the  freak  to  peak  amplitude  of  the  Rayleigh  pulse  of 
the  regional  synthetic  can  he  directly  reflated  to  the  .Ms  valiK'  determiru'd  for 
a  theoretical  Rayleigli  vvave  train  frropagated  out  to  10°  .  The  ri'lat ionship 
between  the  data  fx'ak  to  peak  amplitiidi*  and  its  indirect  .Ms  is: 

Ms(data)  =  .Ms(synth|,oo)  +  log, „((1‘I‘AU„,. )/( IM*.\|s,, 

when'  IM*.\  is  the  peak  to  p<’ak  amplitude  of  the  Rayleigh  pnlsi'.  A  path 
correction  may  Ix’  included  on  the  right  sid<' of  this  expri'ssion. 

This  fialli  correction  is  the  diffen'iice  betwi'r'ii  tlx  '  idividnal  jrath  syn¬ 
thetic  derived  .Ms  and  the  average  theoretical  Ms  for  the  entire  network,  for 
each  .sonr<  <’  reciever  pair,  a  Ms  is  calculated  from  a  synllx'lic  si'ismogram 
firopagated  to  10°  .  Each  such  syntlx'tic  has  tlx'  same  size  source,  so  idi'ally 
one  w'ould  want  each  M,  value  so  measured  to  be  evfual  in  value.  Yet  this  this 
not  so,  for  <'ach  path's  affective  attenuation  at  the  periods  of  inten'st  may 
he  different,  I  he  difference  between  the  mean  network  Ms  and  a  particular 
receiver  .M,  is  tlx-  [lath  correction.  A  negative  path  i orri'ct  ion  value  imjilies 
that  tin-  theoretical  10°  station  .M*  is  larger  than  the  network  average,  'fable 
1  lists  the  network  path  corrections  used. 

fhe  question  arises,  whether  or  not  it  is  valid  to  use  the  average  Earth 
structure  for  a  fiarticular  fiath  to  propagate  a  surfaci'  wave  to  10°  when  tlu’ 
Earth  nxxlel  is  only  meant  to  refh'ct  the  seismic  properties  of  t  he  I'.arth  for  a 
fiath  that  may  only  Ix'  a  small  fraction  of  this  ilistance.  This  is  fiart iciilarlv 
true  of  the  shortest  fialhs  for  which  the  .seismic  waves  traverse  only  wc'st 
ern,  .North  .America,  an  area  of  relatively  high  altennation  compared  to  the 
lontinenial  craton  and  shii’ld  areas.  A  surface  wave  firopagated  10°  along 
a  characteristii  tectonic  North  American  crust  and  matitli'  model  (N  I  S  to 
Dl'd.  for  example)  for  10°  will  be  much  ntore  attenuated  than  a  wave  firofi 
agated  tlx-  same  distance  through  an  average  structure'  from  NFS  to  tlx' 
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raslrni  seaboard  (NTS  to  S(’P,  for  example).  HeiKe  the  ealeulated  for 
the  N  I  S  to  DlKl  structure  would  be  .smaller  than  the  N  I  S  to  St'P  . 

There  are  several  methods  to  correct  for  this  path  dej)eii(lf'iit  effect.  .\s 
(‘xplained  above  one  may  implement  path  corrections  which  acduint  for  the 
thewetical  difference  in  attenuation  between  paths.  .Anotlu'r  means  is  to 
make  a  mixed  path  structure  which  has  the  appropriate  path  structure  from 
the  source  to  the  acliiai  station  distance,  with  the  rest  of  the  path  out  to 
'10°  being  a  gein-ric  .seismic  velocity  and  attenuation  mo»h  l.  I  he  mi.xerl  j)ath 
approximation  for  the  fundamental  Rayliegh  mode  is  easy  to  implement, 
riie  approximation  is  equivalent  to  assuming  that  the  total  horizontal  en¬ 
ergy  flux  is  constant  across  the  transmission  boundary  ( llarkrid<'r.  1981). 
(Levshin,  1985)  and  (Bache  c/  al  ,  1978).  Finite  eleim'nt  results  show  that 
even  for  a  continental-oceanic  crust  transtion  zone  this  approximation  is  rea¬ 
sonable  (llarkrirler,  1981).  For  the  eases  in  this  stinly  where  the  structures 
which  comprise  the  mixed  path  are  both  continental  struc  tures  [i.c.  not  too 
dissimilar)  the  approximation  is  robust  enough  for  the  synthetic  seismogram 
calcualtions. 

We  have  implemented  both  procedures  individually  and  in  c onjutict  ion 
to  se<'  what  thcnr  effects  are.  Another  method  would  be  to  inclinh'  empirical 
station  corrections  (Yacoub,  1983),  ((livcMi  and  .Mc'llman.  198()).  The  findings 
concerning  the  path  corrections  are  discu.ssed  in  the  results  section. 

4  Data  Analysis  and  Results 

riie  .seismograms  were’  band-passed  fillerc'd  between  (i  and  lOO  sc-conds  to 
minimize  contaminating  noi.s«'  as  clescribed  pre\  iously.  I  hi-  \<'rtical  records 
were  visually  ins|)ected  t«)  insure  that  tln-y  were  within  the  correct  time 
winch)W  ancl  that  their  signal  to  nc>ise  ratio  was  above-  2.0.  .\K  values  were- 
then  calculated  for  the  data  as  per  the*  mc-thod  clcscribecl  above  (c'c|natic»n 
10)  with  several  variat  ions.  I  he*  .synt  hc't  i<  sc-ismograms  weie  also  bancl-|)assec| 
filtered  bet  ween  (i  and  100  se< oncls  for  consist c-nev.  I  he-  ,\1,  are  plot  t  ed  against 
seismic  magnitudes  of  scncral  scales  for  the*  same'  se-t  of  evc-nts.  It  should  be 
noted  that  complete  magnitucle  lists  wc-rc*  not  available  for  all  102  e  vents. 

We  chose  to  compare  or  plot  our  data  primarily  with  body  wave'  mag¬ 
nitude’s  de'te'rmiiu'd  by  bilwall  and  McNe’ary  (198')).  Ihc'  l.ilwall  data  sc't 
contains  75  se've'ii  of  the  102  e'vents  exatnine'd  by  us  and  is  belie-vc'd  to  be 


a  well  (IclcrrnirK'd  and  sclf-ronsistani  list  of  fiii,  valiu's  (hat  have  sinall  er¬ 
rors  due  to,  among  oilier  tilings,  the  inclusion  of  network  station  corrections. 
Fig.  2  shows  llu'  nih-yield  relationshi|)  for  ev<-nt.s  in  this  study  for  which  ni|, 
and  yield  information  were  available.  It  is  important  to  notice  that  ev<Mits 
al)ove  and  below  the  water  table  separati^  into  two  distinct  populations.  F'or 
this  data  set  tliis  separtalion  is  only  a|>|>ar<‘nt  near  the  cluster  of  <’vents  with 
tin,  s  aromul  u  l-  .\lso  not  ice  t  lu' v<‘ry  small  error  liars  for  I  his  data;  for  many 
I'vi'iils  the  error  bars  ar<'  smaller  than  the  sMiibols  ilemarking  a  data,  point. 
1  he  corrc'lalion  between  mi>  and  yi<‘l<l  is  good,  with  the  scatter  mostly  being 
due  to  tlu'  above  water  table  shots. 

Several  sets  of  synthetic  Hayleigh-waves  wcne  lah  iilatc'd  at  the  10°  dis¬ 
tance.  OiK'  set  was  propagated  along  the  single  structure  model  (hereafter 
referred  to  as  the  single  path  case)  which  reflects  the  average  Karth  structure 
between  N  IS  and  a  given  station.  Two  sets  of  mi.xed  (latli  synthetic  seis¬ 
mograms  wer<'  also  generated.  For  that  part  of  the  path  beyond  the  actual 
.source  recu'ver  distance,  out  to  10°.  a  generic  eath  structure  was  used.  I'he 
.N  IS-IISSI)  Karth  stnictiirr'  was  chosen  for  this  gc'iuTic  path  section,  a.s  it  is 
a  relatively  simple  structure  which  g«'nerat<'s  slabh'  siirfaci'  wavi's  and  it  is 
roughly  an  inlermediair'  rang**  station  (di.stanc<‘  <  1200  km),  so  t  hat  its  striic 
ture  can  be  considered  to  be  an  “average”  structure'  for  the  lu'twork.  I'Ik' 
difference  between  tlu'se  two  mixed  fiath  earth  structures  is  in  their  spectral 
attenuation  coefficients,  with  ')ft  being  twice  as  large,  at  a  given  frequnecy. 
for  the  mixed  [lath  2  ca,se  as  for  the*  mixed  path  I  case'. 

Siirfae  e'  magnitiiele's  we-re' first  calculated  from  t  he'  10°  synt  hetics  geiie'iate'el 
with  a  single'  structure-  jrropagatiem  path.  Fig.  '.i  elisjilays  single  path  M, 
values,  calculate-d  as  de'scribe-el  abeive,  versus  be)ely wave-  magnitiiele  (iiii,). 
riiese  mi,  s  are-  tho.se  Lilwall  and  McNeary  (1085).  In  the  left  hg.  the-  Ms’s  are 
calciilateel  witheriit  path  corrections,  wlu're'as  path  eorre'ctieuis  are  includf'd  in 
the  figure  on  the  right.  Flie  seilid  line  is  the  best  fitting  weighte-d  le-ast-seiuares 
regression  of  the  data,  with  the  weighting  factor  being  inverse-ly  preeportieinal 
to  indivielual  event  stanelard  elcviations.  'I'he  elashe-el  line-s  re'|)re'sent  the-  twei 
stanelarel  ele-viation  e-rreir  of  the  fit  of  the  line*  to  the-  elata.  Sediel  blae  k  circle's 
are  shots  below  the  wate-r  table,  shots  above  the  wate-r  table'  are  ope-ii  circle's, 
anel  open  sepiare's  are  she>ts  for  which  this  information  is  not  kiiejwn.  Note 
the'  e-rrejr  bars  are-  ap|)re)ximately  50  pe're:enl  large'r  feir  the-  uneorre'cte'el  Ms's 
(Fig.  da)  than  for  the  case  e)f  pat h-correcteel  Ms's  (I'ig.  db).  I'he  scatte-r  in 
the  data  is  alse)  slightly  le-ss  ferr  the  path  corre'cte-el  Ms's.  H  appe-ars  that  the 
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path  corrections  do  improve  M*  measurments. 

The  most  significant  effect  of  including  path  corrections  is  the  reduction  in 
variance  of  individual  magnitudes.  Without  path  corrections  the  imlividual 
station  magnitudes  have  a  hi-modal  distribtition  relied ing  the  t  wo  generic 
Karth  models  of  North  America:  the  tectonic  w«'stern  and  cratonic  eastern 
crust  and  upper  mantle  structures.  The  path  corrections  l)ring-in  the  out¬ 
lying  station  magnitudes  values  towards  the  tnean  value.  I'ahle  1  lists  tlie.se 
network  path  corrections.  'I'he  first  column  lists  the  corrections  for  single 
path  synthetics.  A  positive  value  denotes  that  the  Ms  for  a  station  is  smaller 
than  the  network  theoretical  average.  The  fourth  column  is  the  number  of 
events  that  W('r<'  recorded  at  the  station  for  tin*  flata  ('iitire  ,s<'t . 

W('  next  explored  t  he  <'ffect  of  mix«'d  path  transler  fuiu  tions  ii|)on  the  .Ms 
calc.ualtions.  As  described  above,  we  chose  the  path  to  HSSl)  as  a  generic 
structure  for  the  second  portion  of  the  mixed  path  synthetic  seismogram 
calcuations.  We  generated  two  sets  of  these  .synt lu'tics.  I'Ik'  attenuation  of 
the  generic  |»ath  was  doubled  for  one  of  tlu-se  .sets  (lfS.Sl)x2).  !•  ig.  1  shows 
th<'  attenuation  factors  (gamma)  as  a  function  of  period.  I  lie  liiu'  labeled 
RSSDx'i  is  that  of  the  increased  att<‘ntuation  stnuture.  It  is  referred  to 
as  “mix<-d  path  2”  througliont  this  stinly.  I’lu'  l(»w('r.  dashed  ( iirve  is  tlu' 
attenation  curv<'  for  tlu*  HSSl)  struct nr<’.  SyntIu'tics  made  with  this  HSSl) 
generic  structure  for  tiu'  latter  portion  t)f  the  10°  t raved  path  will  he  re  hure'd 
to  iis  “mixed  path  1”.  Fable  I  gives  the'  path  corree  tions  IVir  eac  h  station  foi 
these  two  cases,  also. 

Fig.  5  is  analogous  to  Fig.  3,  tlie  difference  being  that  the  .\1^  magnitudes 
were  calculated  using  synthetic  .seismograms  using  tlu'  mi.\('<l  path  1  :nofl(d. 
In  Fig.  ."ja  the  Mj's  are'  calcualte'd  without  path  eeirre-ct iem  te'rms.  while'  in 
Fig.  .51)  path  ce)rre'e'tie)ns  are-  incluele'el.  .\s  befeue-  the'  addition  of  the'  path 
correction  terms  cut  the'  varianee'  by  aleeuit  a  fae  te)r  e>f  2.  .Mso  tin-  a\e'rage  .M, 
value' ele'eTe'ase'el  by  a  0.01  units  freun  lhe>se'  witheuit  path  eor  rce  I  ions,  lids 
ehlfe'ce’iiee'  may  be'  st  at  ist  ie  all  v  sigidlieant ,  le>r  t  h<'  si  a  in  lard  on  or  oi  I  lu'  me 'an 
e»l  the-  re'gre'ssiem  line-  is  einlv  0.02(i.  What  is  most  sliikiinj  hclwi-on  hi';  s  t 
auel  .h  is  that  slojee'  eif  I  he’  re'gre'ssieui  line'  are  t  in'  same’  her  tin'  i  \ve>  mixeei  pat  h 
I  ease’s,  with 

mi,  =  0.8:1  X  M,  d-  M. 

For  the’  single’  pat  h  e  ase’,  pat  h  eeure’cte’el  M„  value’s  give’ t  he’  same  i  elat  ionshij). 
but  the’  slerpe  is  appre’eiable’v  large’r  (  0.0  )  feir  the’  mie  eure’e  le’el  magnitudes. 
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although  I  ho  (lifToronco  lios  within  tho  orrors  l)oiin<is.  It  would  sooin  that 
both  path  rorroct ions  aiul  and  inixod  path  (Inx'u's  funictions  iniprovr  Ms 
d(‘t<Tminat ions  for  the  nn'thod  us«‘d  horo. 

I'hus  th«‘  inclusion  of  path  corr<'ct ions  in  <al<  ulaling  Ms  1ow(ts  tlu*  mag¬ 
nitude  value  by  0.01  units  on  averag«-.  I’his  can  be  (‘xplained  by  the  network 
( ()V('rag«'  atid  t  he  range  of  syntlu‘ti<-  .M,  valiH's.  The  dist  ribiil  ion  of  Ms  values 
is  skewe<i.  with  th<'r<‘  Ix'ing  a  signilicant  number  of  stations  with  llu'oretical 
Ms  values  signiiicantly  larger  than  most  of  the  stations  within  tlu*  network. 
Kvents  for  which  such  stations  reported  would  yield  a  largi'r  average  M*  than 
for  events  1  hat  did  not . 

Fig.  (ia  and  (ilj  ar<'  .Ms  vs.  nii,  plots  for  tlu'  mixe'd  path  2  case  without 
and  with  path  c(jrr('ctions.  respectivxdy.  For  (ia  a  fixc'd  slope  regression,  was 
use<l.  so  that  more  a  more  ai)t  comparison  could  1h'  made  to  its  mixed  path 
1  (ounterparl  |)lot.  Comparing  Fig.  to  Fig.  (ia  shows  (hat  the  average  Ms 
value  droi)|)ed  0.02  magnitude  units  for  the  mixed  path  2  cas('  r<'lativ(>  to  the 
mixed  path  1  case.  This  follows  from  the  fact  that  the  attenuation  for  tin* 
second  ()or(ion  of  the  path  is  twice  as  large  for  (he  mix('d  path  2  rase  as  that 
for  the  mixed  path  1  case. 

In  comparing  Fig.s  5b  and  Gb,  however,  it  is  a|)|)arent  that  there  are  no 
rlilferences  in  magnitude  valu<-s  and  the  regression  line  in(('rc('pt  and  slopes 
ha ve  I  he  same  valu<- in  both  cases,  even  though  tin*  mix<'fl-pal  h  (Ir<'en’s  func¬ 
tions  used  f<jr  .Ms  cahulalion  purposes  dilhued  wuth  resfx'ct  to  attc'imation 
for  the  two  (  ases.  I  he  effect  of  the  path  correction,  besides  reducing  the  data 
variance  as  described  above,  is  to  negate  the  effect  of  differences  in  attenua¬ 
tion  betwfen  the  tw'o  models.  To  obtain  stable,  robust  M,  values  with  this 
method  it  is  best  then  to  use  mixed  path  generated  synthetics  in  conjunction 
with  path  corrections  for  the  10°  Ms  measurements.  I'he  variance'  among  the 
mixed  path  based  M*  values  for  the  network  is  smaller  than  that  when  Ms  is 
derived  frmn  single  path  syntln’tics,  so  (hat  magnitude  nK'asurenu’nts  will  be 
more  consistent  when  tlx’y  are  determined  from  mixed  path  synthetics.  I'liis 
is  particularly  imi)ortant  for  events  with  few  reporting  stations.  All  further 
plots  of  .Ms  in  this  study  use  values  obtained  from  the  mixed  path  1  case  with 
path  corrections,  unless  stated  otherwise. 

How  w<'ll  the  final  Ms  values  reflect  the  actual  seismic  magnitiuh'  of  these 
events  ne< fssitates  having  another  measure  of  tlu'ir  si/.('.  In  the  <'vent  of 
anomalously  high  or  low  seismic  source  coupling,  for  example,  both  body 
waves  and  surface  waves  should  be  affected  similarly  by  tin'  coupling  effect. 
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A  magiiilud*'  |)araiii<‘lcr  iiidrpoiit  of  srismir  of  oliscrvat ions  would  ho  iiso 
fid  to  plot  IIk’  agfiinst,  so  wo  liavi'  also  (illcd  oiir  rosiills  lu  oslinialod 
log  yiolds.  Fig.  7  shows  tin*  rolationsliip.  Yiold  valuos  arc  ostimat<‘d  to  ho 
within  10  percent  of  the  actual  yield  (Springer  and  Kinnaman.  11)71).  ic'ld 
information  wa.s  available  for  97  of  the  events,  thus  yields  mak('-u|)  the  most 
comprehensive  data  sot  to  compare  our  results  to  as  well.  I'lio  yiolds  for  this 
data  set  range  over  throe  onlers  of  magnitude  in  size.  I'ho  grt'atost  scatter, 
as  in  the  case  of  Mb  vs.  log  yield,  is  due  to  shots  above  the  water  table.  It 
should  also  be  be  kept  in  mind  that  the  scattr-r  would  bo  further  reduced 
if  tlw'  data  were  separati'd  into  populations  basivl  on  their  location  at  .NTS. 
(j.c.  Pahute  Mesa,  Rainier  Mesa  and  Yucca  F-lat).  Bc'cause  of  the'  classified 
nature  of  some  the  yields,  it  is  not  possible,  here,  to  closely  e.xamine  the.se 
efh'cts  with  respect  to  yield. 

Since  our  magnitude  values  are  ba.sed  on  thc-orelical  continental  struc¬ 
tures,  as  well  as  tin;  particular  network  used,  we  wanted  to  compare  our  M* 
values  to  those  obtaiiu'd  from  standard  M*  methods.  One  such  standard  data 
s<'t  is  that  of  Marshall  cl  al  (1979).  'I'here  is  a  overlaj)  of  11  events  Ix'twt'ii 
studies.  We  preformed  a  fixed-slope  (slope=1.0).  rc'gression  of  our  Ms  values 
to  theirs.  Fig.  8  shows  that  the  correlation  is  viuy  good:  siattc'r  is  small  for 
I'vents  above  and  below  the  water  table.  It's  important  to  note  that  with 
our  method  we  are  able  to  tiieasure  M*  for  evi'iits  0.75  units  smaller  than  the 
smallest  Marshall  values.  \V<-  are  able  to  measure  .Ms  for  these  smaller  evc'iils. 
because  we  are  able  to  make  use  of  near-regional  (  <  500  km)  records  with 
tlu'  m<'t  ho<l  described  in  this  paper.  I'Ik' offset  in  .M,  bet  weiMi  si  ales  is  0.5.'f 
with  mean  standard  error  of  O.Od  magnitude  units.  This  olfset  is  due  in  part 
to  the  difference  in  definition  of  Mj.  .At  10°  the  offset  in  magnitude  is  0.15. 
thus  reducing  the  olfset  to  O.dS  units.  However,  for  the  method  decribi'd 
in  this  study,  Ms  is  baseil  u|)on  a  theoretical  ni'lwork  average  .'.F.  so  it  will 
have  a  bias  attached  to  it  which  is  de(>endent  u|)on  the  netwoik  usi'd.  This 
netw’ork  bias  can  be  a.ssumed  to  be  ies|)onsible  fur  jiart  of  the  olfsi't,  as  well. 

I’abh’  2  lists  th<’  linal  mixed  palh,  juith  corrected  .M^  values  for  the  102 
events  of  thisstudv.  1  he  first  column  lists  the  nnmliei  of  si  al  ions  recording 
the  I’vent.  NexI  ari’  given  I  he  surface  wave  magmiude  and  associated  erim 
for  the  event  are  given.  Next  is  a  three  letter  shot  informalioii  i ode.  I  he 
first  letter  denotes  its  geographic  loialion:  5  ucca  (Y).  Pahnie  (P).  Rainer 
(R),  or  Climax  Stock  (C).  I  lu’  s<‘cond  is  whether  its  shot  de|)tli  was  above 
(A)  or  Ix'low  (H)  the  wafer  table.  'I'lu’  last  letter  desi  ribes  the  shot  site  rock 
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as  tuir  (  1  ),  rhyolite  (1{).  granite  (C),  or  alhivinm  (A).  An  uiulesrseore  means 
that  th<'  information  is  not  known.  Tlie  final  two  roliimiis  are  tlu^  event’s 
name  and  Julian  data,  respectively.  The  events  are  listed  in  chronological 
order. 

To  determine  the  portahiltiy  of  this  Mg  calcnlation  method  tin*  events 
need  to  h<'  se|)arated  into  groups  based  on  their  source  regions  and  then 
compared,  one  group  to  another,  in  order  to  se«'  if  there  are  .systematic  dif¬ 
ferences  in  M„  valiK's  relativ<'  to  any  other  magniliuh*  scale.  'I'hrer'  main 
g(’ographic  source  regions  <<)mpris<‘  the  «‘vent  <lata  set:  I’ahule  Mc'sa,  Kainier 
.\I('sa  and  't’licc  a  I'  lat . 

Whether  or  not  a  shot  occurs  within  saturated  material  is  anotlu'r  crite¬ 
rion  by  which  to  separate  events  in  order  to  look  for  systcunatic  differences 
in  .Ms  valiK's.  Other  studies  have  found  significant  seismic  coupling  differ¬ 
ences  betwisMi  (’.xplosions  detotiated  above  and  below  tlu-  water  table  (Cupl  a, 
1989).  .so  it  a  reasorial)le  parameter  to  study.  Reviewing  Fig.s  2  and  7,  it  is 
also  a|)parent  that  for  shots  fired-off  below  the  water  table  have  a  larger 
seismic  magnitude  than  those  <ietonated  above  the  water  tabh'. 

I'ig.  9a  shows  (he  relationship  between  M,  vs.  Lilwall  \\\),  for  all  N'l'S 
events.  The  siirfac('-wave  maginitudes  were  all  calculated  using  mixed-path 
Green's  functions  and  i)ath  corrections.  Fig.s  9b  and  9c  divide  the  data 
[><;pnlatioris  into  above  and  below  the  water  table,  respc'ctively;  shots  for 
which  water  table  information  was  not  available  w('re  left  out.  'I'here  is  no 
appreciable  difference  between  the  above  water  table  and  below  water  table 
curves.  I  his  is  not  surprising;  refering  to  Fig.s  2  and  7  it  is  appareiH.  that  bot  h 
seismic  magnit  inh'-yield  (  urv(*s  show  that  for  a  givuui  yield  an  event  below 
the  water  tabh’  has  a  larger  magnitude  than  a  shot  above  the  wat<'r  table. 
It  follows  that  the  .Mg  -  nib  relationship  may  not  show  the  same  discrepancy 
betwe<-n  shots  detonated  above  and  below  the  water  table,  for  the  effect  of 
the  water  content  in  tin'  shot  medium  shouhl  affec  t  surface  waves  and  body 
wavesin  the  sann'  manner.  .All  three  regression  curves  are  es.sentially  the  same 
within  t  he  error  bounds.  There  is  considerable  scatter  in  all  three  figures,  but 
that  is  not  surprising  considering  the  diversity  of  the  samjiled  populations. 
Fven  with  this  scatter,  the  best-fitting  Mg-rnb  liin*  is  well  constrained,  for  the 
population  covers  a  wide  range  of  magnitudes. 

Fig.  lOa  gives  the  Mg-inb  relationship  for  all  ’i'ucca  events.  I'he  rc'grc'ssion 
curve  is  signifcantly  different  from  that  of  Fig.  9a.  I  he  scatter  in  the  data 
is  reducc'd  by  in  percent  over  that  of  the  general  population.  Sejiarating 
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Ulo  events  witli  respect  tlieir  relation  to  the  water  tal)le  yields  two  distinct 
curves,  unlike  the  case  where  shots  from  all  sites  are  grouped  together,  lor 
the  Yucca  events  below  the  water  table  (Fig.  lOb)  the  regression  curve  lit 
is,  within  the  errors,  is  not  greatly  different  from  the  case  (jf  all  .\  l  .S  events 
below  the  water  table  (Fig.  9b).  In  Fig.  9c,  for  shots  in  dry  nudium,  the 
Mg-mb  curve  is  significantly  different  from  the  NTS  above  water  table  curve 
(Fig.  9c).  The  significance  of  this  curve  is  questionable,  however,  due  to  the 
paucity  of  data  used  to  establish  it. 

h'ig.  11a  plots  all  Pahute  event  Mg’s  vs.  their  respective  nib's.  Th(>  re- 
altionship  is  essentially  the  same  as  for  the  aggregate  N  TS  plot  (Fig.  9a). 
Fig.  Ill)  shows  the  relationsliip  for  Pahute  shots  below  tin'  water  table.  1  he 
slope  of  this  curve  is  nearly  the  same  as  that  for  Yucca  Ixdow  water  lalile 
shots,  although  the  intercept  differs  appreciably.  This  result  implies  that  for 
a  given  mb,  surface  wave  magnitudes  for  events  at  Pahute'  Mesa  are*  larger 
than  those  at  Yucca  Flat.  Fig.  11c  shows  the  Ms-nii,  relationshii)  for  Pahute 
events  above  the  water  table.  Because  events  in  this  category  are  clustered 
around  rrii,  =  5.5,  a  well  constrained  line  cannot  be  obtained,  .so  we  applied 
a  fixed-slope  regression,  using  a  slope  of  0.82  (that  of  tin'  below  water  ta¬ 
ble  case),  (’ompariiig  Fig.s  lib  and  11c,  an  off-set  in  Ms-mi,  of  0.09  with  a 
mean  standard  error  of  0.04  is  obtained,  implying  some  degree  of  difference 
in  seismic  cou[)ling  between  body  and  surface  wave's. 

Fig.  12a  displays  the  M,-mb  regresssion  curve' and  data  for  llainie'r  .Vle'sa 
events,  all  of  which  are  above  the  water  table.  Although  the'  <  lustering  of 
data  near  nii,  =  5.0  cause's  the  curve  te)  be  poorly  ce)ustraine'd,  a  slope  cleese 
te)  that,  for  Pahute'  e've'iits  anel  Yucea  (b<'le)w  the'  wate'r  table')  is  e)btaine'el. 
(Mnstraining  t  he  sle)p<'  to  be  0.82,  as  in  the'  e  ase'  of  Pahute  e've'nts  yie'lds  1' ig. 
121).  (om[)arisie)n  eif  the-  <'e|uatie)n.s  e)f  Fig.s  12b  and  11b  give'  an  olf-se't  e)f 
O.ill,  with  a  mean  standarel  e'rre)r  e)f  0.04,  nib  units,  so  that  for  a  given  nu,. 
slmts  at  Pahute  proeluce  larger  Mg  value's  than  at  Hainie'r  .Me'sa.  Fithe'r  the 
Pahute  site  is  meue  e'flie  ie'iit  at  preielucing  surface  wave's  or  the'  Hainie'r  site'  is 
le'.ss  e'llicieut  at  eoupling  beidy-wave  ene'rgy.  Hainie'i  e've'iits  are  tunne  l  shots, 
riie  immeeliate'  seuirce'  re'gieui  (H  <  200  m)  may  be'have'  like'  an  asymme  trie 
cavity,  re'sulting  in  a  .seiurce  that  is  nein-iseitreipie-  (Zhao  and  llarkride-r.  1991 ) 
anel/or  .seismic  eeiupling  that  has  streing  freque'ticy  de'|)e'nelene e. 
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Conclusion 


1  lie  meth(j(i  df'scrihed  li<'r<'iii  to  <al(iilat<'  surface*  wav'c  tiiagtiitudes  allows 
tlie  ineasniing  of  for  micl<*ar  <*x|)losioMS  ove'r  a  wide*  iiiagiiil  iid<*  <lis(.ril)ii- 
lioti.  I.sing  tills  lecliiii<|Ue,  it  is  now  possilrU*  to  use*  near  n'gioiial  (A  <  ITi®), 
long  period  records  in  con jiiiu  t ion  with  tin*  far  regional  (A  <  15°)  and  I.ele 
seisiiiie'  re<  ords,  that  have  previoiislv  l)<‘<“ti  us<*d  to  make'  (  (invent ional 
ineasiin'iiK'nl s.  in  order  to  measure*  snrfaee  wave  magnit nele*s.  I’liis  iiirre*ase* 
in  ohservatieins  has  two  advantages.  First,  for  any  e've'iit  the  statieui  netweirk 
coverage  is  enhanecel  in  terms  eif  ove*rall  numl>e*rs  as  we'll  as  in  a/.imuthal 
coverage,  in  particular  stations  only  several  hundred  miles  away  from  N'l'S 
in  ('alifornia.  .\e\ada  and  I’tah  can  be  include*el  in  a  network  that  otherwise* 
would  have  nei  ceiverage*  to  lire  w(*st  or  southwest.  rh(*se  improve'nu'uts  make 
the  network  .\h  s  more*  stable*  and  statistically  robust.  Se*conelly.  smaller 
e*ve*nts  with  surface*  wave*s  that  haven’t  been  analy/('el  will  now  have  such 
observations  available,  ser  that  the*ir  surface  wave*  magriituele*s  now  can  be* 
calculated,  I’lie  m(*thoel  proeluees  stable  M*  value's  that  are*  consistent  with 
other  seismic  magnituele*  scales.  Thirdly,  the  eirect  erf  inaccuracu*s  in  estimat¬ 
ing  Q  are'  rie*gligible  for  very  near-re*gional  recordings  (<  500  km).  With  the 
M,  calciil.il  ioii  t('chni(|U(’  use'd  lu're*,  one*  can  take*  advantage*  of  sue  h  ne*arby 
recordings,  l.astly  this  iii(*t)iod  make*s  it  ve*ry  e*asy  to  use  liislorical  analog 
data  s(*ts  iiioie'  e'asily.  for  it  is  not  ne'e'e*ssary  to  use*  digit i/.e'd  data  if  only  tlie* 
ma.vimiiiii  pe'ak  to  jie'ak  amplitude  nee*ele'ei  to  ealciilate'  .M*  s  in  this  fashion. 

In  oreler  to  elo  a  thorough  inv(*sl igatiern  of  serurce*  elfe'cts  discus.se*el  above, 
a  more*  ce)mpl(*le*  elata  se*t  of  e*V('nts  with  elijferenl  pe'itinent  souree*  param- 
e'ters  is  n('e*(le'(l  feir  observatiein.  Spe*cilically.  eve'iits  cherse'u  ferr  any  one  site 
lercation  should  cover  a  wide*  range  erf  yields  err  magnitude’s  and  eomprise 
eve’iits  eleternate’d  both  aberveand  below  the  water  table.  Obviously,  at  some* 
site  lercations.  it  may  iiert  be  perssible  ter  come  up  with  such  a  eermprehen- 
sive.  iele*al  data  set.  This  is  jrarticularly  true  erf  Hainie*r  Me’sa  where  merst  all 
(’vents  have’  yie*lds  le*ss  than  'iOKt.  Frerm  the*  re’sults  obtaiiu'd  with  the  data 
set  iiseel  here*,  tlu’re  (loe*s  apfre’ar  ter  be*  signicant  differe'iice's  in  .s(*ismic  con 
pling  be'tW('e*n  .N  TS  sub-sites,  wdth  <*ve*nts  at  I’ahute'  Me'sa  piodiicing  large'i 
sill  face- wave*  magnitude's  ferr  a  give*n  nii,  than  at  Hainie’r  Mesa  err  ^’iice  a  Flat. 

I  his  (liscr(’[iancy  is  largest  ferr  Kainier  Me*sa  eve'uts. 

.\s  stated  e'arlie’r  the*  me*therel  use*el  he*re  ter  obtain  .Ms  value's  le’iids  itse'lf 
to  such  a  diverse’  data  se’t,  ferr  the*  measure*ment  of  the’  smalle'st  e've'iits  is 
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facilitated  with  near-regional  observations. 
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Table  1;  Network  Path  Corrections 


Path 

Correction 
Single  Path 

Path 

Correction 
Mixed  Path  1 

Path 

Correction 
Mixed  Path  2 

Station 

Name 

0.12 

-0.21 

-0.20 

AAM 

0.13 

0.20 

0.20 

ALO 

0.52 

0.49 

0.49 

BKS 

0.19 

-0.40 

-0.40 

BLA 

-0.36 

-0.12 

-0.20 

BLC 

0.29 

-0.02 

-0.02 

B02 

-0.47 

-0.30 

-0.29 

CMC 

-0.17 

0.03 

0.03 

COR 

0.08 

-0.39 

-0.40 

DAL 

-0.08 

-0.07 

-0.07 

DUG 

0.29 

-0.10 

0.06 

EDM 

-0.04 

-0.19 

-0.19 

EPT 

-0.27 

0.02 

0.02 

FCC 

-0.27 

-0.06 

-0.07 

FFC 

0,19 

-0.26 

-0.26 

FLO 

-0.36 

-0.31 

0.20 

FRB 

0.40 

0.42 

0.41 

FSJ 

0.19 

-0.25 

-0.26 

FVM 

0.19 

-0.24 

-0.23 

GEO 

-0.15 

0.01 

0.01 

GOL 

0.24 

-0.09 

-0.09 

GSC 

-0.44 

-0.19 

-0.19 

INK 

0.09 

0.16 

0.16 

JCT 

-0.18 

-0.06 

-0.06 

LHC 

-0.12 

-0.07 

-0.06 

LON 

0.33 

-0.21 

-0.01 

LUB 

-0.47 

-0.42 

-0.44 

MBC 

0.12 

0.31 

0.30 

MNT 

0.29 

-0.04 

-0.03 
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-0.28 

0.10 

0.16 

OGD 

0.12 

0.34 

0.34 
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-0.10 
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OXF 
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0.43 
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0.15 
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0.14 

-0.08 
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-0.02 
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STJ 
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0.14 
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0.19 
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-0.34 

YKC 

0.13 
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RSNY 

-0  44 
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0. 12 
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GAC 
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No-  Ms  Shot  Julian 

-'---i-  Ms  s.d.  Info.  Event  Date 
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Figure  1..  Map  of  North  American  station  network  useii  in  tliis  study.  Tlie 
"spoked  wheel"  is  the  .Nevada  Test  .Site. 
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mb(LW)  vs  Log  Yield 


Log  Yield 

(s.d.=0.16) 

Figure  2.  Lilwall-iiii,  vs.  log— Yield  for  events  from  this  study.  Th'  '-wli  1  liu. 
the  best-fitting  regressit)n  line.  The  d:ished  lines  show  tin-  j  « nent  r..nri  i.  ii  c 
interval  of  this  line.  Lines  through  the  data  points  rei)res<  nt  one  standard  devia¬ 
tion  in  a  datum  mea-surement.  Blackened  cirlces  represent  sourcis  iKiuath  tin 
water  table,  open  circles  arc  events  above  the  water  table  and  ojicn  sf)uares  are 
events  for  which  this  information  is  not  known.  Events  detonated  below  the  wa¬ 
ter  table  have  a  larger  111^  for  a  given  yield.  Besides  this  separation  of  data,  tli'  i" 
is  little  scatter  to  the  data.  The  consistency  of  the  ni( -yield  relationship  make-  ii 
reasonable  to  use  these  nij,  values  to  plot  our  M^,  rnea-surements  against. 
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Attenuation  vs  Period 
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Figure  4.  'I'lic  l\><)  :il tcinKilioii  tiiodcis  lur  tlic  gciirric  |)(jrii'»ii  ol'  the  li\ liiirl 
pagalion  modt‘1  arc  plutlcd  v('ri>us  period.  Gamma  is  the  atteuuatiun  eoeUKlrm  ni 
a  given  period.  Model  KSSDx'J's  attenuation  is  twice  tlial  of  tie  HSSl)  mol'!. 
Mj,  values  ealeiilat<'(l  with  synthetics  using  model  K.^Sl)  are  referred  to  as  mix'  I 
path  I,  while  values  determined  from  synthelies  created  using  attenuation  mode] 
|{.SSDx2  are  referred  to  its  mixed  path  2. 
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Figure  7.  vs.  log  y'ldd  for  all  available  NTS  evenst.  Surface  wave  niagiii- 
tudfs  were  calculated  with  mixed  path  1  Green’s  functions  and  path  corrections. 
Scatter  is  mostly  due  to  differences  in  shot  medium  and  source  region. 
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Ms(WH) 


Ms(WH)  vs  MS(HSR) 


Ms(MSR) 

Ms(WH)  =1.00  X  Ms(MSR)  +  0.57 
(s.d.=0. 11) 


Figure  8.  A  conipari.soii  of  our  values  to  those  of  Mar  hall  d  nl.  (I'.oO)  h  r 
I  i  events.  A  slope  of  1.0  was  assumed  for  this  regression.  Mth'  imh  there  is  an 
off-set  of  O.rj.'i  units,  the  relationship  is  very  consistent. 
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the  water  table  and  (c)  events  above  the  water  table.  The  best-fittiii!;  linear  re¬ 
gression  model  of  the  data  is  plotted  a'^  a  solid  line,  with  it's  eciuation  printed 
below.  The  dashed  lines  are  the  two  standard  deviation  confidence  intervals  of 
the  regression  line. 
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Figure  11.  Mj  vs.  Lilwall  m^,  for  (a)  all  available  PaluUe  Mrsa  events,  (b)  events 
below  the  water  table  and  (c)  events  above  the  water  table.  ’I'lie  best-Htting  linear 
regression  mode!  of  the  data  is  plotted  a-s  a  solid  line,  with  it’s  equation  prime. | 
below.  The  dashed  lines  are  the  two  standard  deviation  coiilidence  intervals  of 
the  regression  line. 
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A  Note  on  Hie  Relocation  Tibetan  Earthquakes 

Lian-Sht  Zhao  and  Donald  V.  Htlmberger 
Seismological  Laboratory,  252  -  21 
California  Institute  of  Technology 
Pasadena,  California  01125 

Abstract  Present  ISC  locations  for  earthquakes  beneath  Tibet  indicate  a  random  distri¬ 
bution  of  events  down  to  a  depth  of  about  50  km.  This  distribution  might  be  expected 
from  a  relatively  cool  crust  which  would  allow  the  seismo-genic  zone  to  extend  to 
such  depths.  A  detailed  investigation  of  the  Tibetan  earthquakes,  with  magnitude 
greater  than  5.5  from  1964  to  1986,  yields  a  distinctly  different  picture.  Waveform 
modeling  of  depth  phases  such  as  pF  indicates  that  only  three  or  four  events  from  this 
population  is  actually  deeper  than  25  km.  These  few  events  occur  near  the  edges  of 
the  Plateau  where  active  subduction  is  occurring  as  suggested  by  the  thrust-like  nature 
of  their  mechanisms.  The  events,  averaging  the  entire  population,  occurred  earlier  than 
indicated  by  the  ISC  by  about  3  seconds  which  leads  to  about  a  1 .5%  and  0.59r  over 
estimation  of  F„  and  5„  velocities  respectively  applying  ISC  tables  and  standard  flat- 
layered  models.  A  more  serious  error  occurs  if  the  F„  and  5„  velocities  are  deter¬ 
mined  by  correcting  for  source  depth  but  assuming  the  ISC  origin  times. 

Data,  Analysis  and  Results 

The  events  studied  are  listed  in  Table  1.  Included  in  this  table  are  the  ISC 
determinations  of  depth,  location  and  origin  times  of  these  events  along  with  new  esti¬ 
mates  obtained  from  this  investigation.  The  main  objective  of  this  report  is  to  docu¬ 
ment  the^e  results  and  discuss  their  significance. 
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The  data  comes  from  the  WWSSN  chip  library  of  the  SeismologicaJ  Laboratory 
at  California  Institute  of  Technology  and  consists  mostly  of  long  period  P-waveforms 
and  travel  time  picks  read  from  the  short  period  vertical  components.  The  P-waveform 
riata  is  used  to  constrain  depths  while  the  time  picks  are  used  to  relocate  and  determine 
the  origin  times. 

Errors  in  depth  determinations  using  travel  times  become  serious  when  no  local 
stations  are  available,  which  is  the  case  in  Tibet  Essentially,  there  is  a  complete 
trade-off  between  origin  time  and  depth.  This  difiBculty  can  be  circumvented  by  deter¬ 
mining  the  depths  independently  by  applying  waveform  modeling  of  the  P-waveforms 
{  Langston  and  Helmberger  1975  ).  This  technique  models  the  interference  between 
direct  P  and  free  surface  reflections  to  fix  the  epicentral  depth  although  it  is  dependent 
on  the  velocity’  model  of  the  source  region.  For  example,  assuming  a  surface  layer  of 
velocity  5  km^ec  instead  of  5.5  km^ec  introduces  a  -1.0  km  depth  error,  if  the  travel 
time  (Sfference  between  P  and  pP  is  used  in  modeling  an  event  at  a  depth  of  10  km. 
Thus,  the  source  depth  as  determined  by  P  waveform  modeling  is  not  very  sensitive  to 
crustal  velocity  and  can  be  used  to  estinaate  the  depths  of  the  above  events. 

Many  of  the  P-waveforms  from  Tibetan  events  have  been  modeled  previously 
by  a  number  of  authors  as  indicated  in  the  Table.  These  solutions  have  been  checked 
by  adding  station  coverage  and,  in  general,  found  to  be  satisfactory.  New  depth  esti¬ 
mates  from  modeling  these  events  assuming  the  crustal  model  given  in  Zhao,  Helrn- 
berger  and  Harkrider  (  1991  ),  model  TIP,  are  listed  in  Table  1. 

We  used  the  P  wave  travel  times  for  relocations.  These  usually  agreed  with  the 
ISC  reports  except  in  a  few  cases.  More  than  30  travel  time  picks  distributed  in  az¬ 
imuth  are  available  for  most  events,  thus  providing  excellent  coverage.  The  station 
coverage  is  indicated  in  the  Table  1.  The  relocation  procedure  used  is  a  simple  itera¬ 
tion  scheme  based  on  the  definition  of  the  ray  parameter  p ,  namely. 


P* 


dR 


or 


d/f«= 


P 


(1) 
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where  dT  is  a  small  chaxige  iii  trav  el  time  associated  wiLj  a  smaJl  cliong'^  iii  iocaiioii. 
We  start  with  the  origin  time  and  location  given  by  ISC.  We  define 


dT~{Ta.-To)-T,  (2) 

where  is  the  observed  time  pick,  7^  is  the  origin  time  and  7^  is  the  theoretical 
travel  time  based  on  an  earth  model.  We  assume  TIP  is  the  appropriate  model  at 
source  region  along  with  elevation  correction  and  the  model  JB  at  the  receiver  stations. 
Station  corrections  were  applied  from  the  Dziewonski  and  Anderson’s  (  1983  )  station 
correction  table.  By  assuming  that  the  dTs  are  from  a  Gaussian  distribution,  we  form 


the  expectation  estimate 


1 

6T-^  —  '£dTi. 


and  obtain  a  new  origin  time  Tonew^TQgi^+ST.  From  equations  (2)  and  (3),  we  can 
see  that  5  T  based  on  this  new  origin  time  is  zero. 

Then  we  use  the  new  dT  assuming  TQ^tw  given  in  Equation  (2)  to  determine 
the  nuslocation  dR  from  equation  (l).  We  find  the  estimated  mislocalions  in  the  North 
dx  and  East  directions  dy ,  namely, 


dr* - ^dR^cos^^ 

”i-0 


</y  = - Vd7?,sin4>, 

”  i-C 

assuming  Gaussian  distributions  of  the  two  direction  mislocalions.  where  is  L*.: 
azimuth  of  the  source  to  the  station  (i).  The  new  location  of  the  earthquake  is  given  by 


latitudes  0^,^  and  longitude  <f>ntu,‘ 


^new*“^oid’^dx/ R£  (5) 

^new 

Next  we  use  the  new  origin  time  and  new  location  as  initial  values,  we  use 
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equations  (3)  and  (5)  to  calculate  the  new  origin  time  and  location,  until  some  criteria 
b  met  The  criteria  we  used  here  is  that  the  difference  of  the  standard  error 


and  the  standard  error  of  the  last  iteration  is  smaller  than  0.01 

".-1 

seconds.  In  thb  inversion  we  assume  that,  the  source  depth  is  known,  but  this  assump¬ 
tion  b  not  required  by  the  method.  The  method  converges  fast  normally  imder  5  itera¬ 
tions  and  the  final  results  is  not  dependent  on  the  initial  input.  Since  changing  the  ori¬ 
gin  time  trade-offs  with  Tg,  the  ISC  locations  are  not  effected  much  by  changing  the 
source  depths. 

After  relocation,  if  the  travel  time  residual  of  certain  station  is  greater  than 
three  seconds,  we  consult  the  record  again.  If  there  is  no  problem  with  the  pick,  we 
use  it;  if  the  pick  is  not  certain,  buried  in  the  noise,  we  may  throw  it  away  depending 
on  the  confidence  of  the  pick,  and  event  size,  and  all  other  station  conditions.  We 
throw  away  data  that  have  a  residual  more  than  5  seconds,  even  though  the  on-set  is 
sharp.  For  example,  UME  and  SEO  produce  more  than  10  second  residuals  for  some 
of  the  large  earthquakes,  and  perhaps  were  caused  by  improper  clock  corrections  or 
reporting  errors. 

We  summerize  our  results  in  Table  1.  In  the  table,  “  reic  ”  means  relocation,  “  D 
"  b  the  distance  between  the  relocated  location  and  ISC  location,  “AT”  is  the  relocated 
origin  time  difference  from  that  of  ISC.  We  also  give  number  of  stations  used  under" 
N”,  Under  “  Data  ”,  we  give  the  references.  Hie  standard  error  of  the  origin  time  due 
to  the  inversion  are  also  given  in  Table  1. 

Figure  1  displays  the  contrast  between  the  ISC  depth  determinations  versus  the 
new  estimates.  The  results  are  dramatic  where  all  the  events  occurring  beneath  the 
interior  of  Tibet  become  shallow  events  while  the  evente  along  the  frontal  thrust  zone 
become  deeper.  This  implies  a  warm  lower  crust  beneath  Tibet  similar  to  other  tec¬ 
tonic  regions  of  the  world.  In  Figure  1,  more  than  90  per  cent  of  the  evente  occurred 
in  the  Tibetan  Plateau,  with  magnitude  greater  than  5.5  since  1964,  are  included. 
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The  mislocations  are  not  very  significant,  the  greatest  is  21.1  km.  The  origiri 


tin.e  toi  ..(iona  are  more  serous.  Of  the  59  earthquakes,  .54  earthquake  ;  '  v-r.  -  ■  > 

Uer  than  reported  by  ISC.  The  average  of  all  earthquakes  is  -3.1  seconds,  or  these 
earthquakes  have  occurred  3.1  seconds  earlier  than  reported  by  ISC.  The  worst  case  is 
-7.4  seconds  for  the  second  March  6,  1966  earthquake.  Eleven  of  the  events  have  a 
residual  greater  than  5.5  seconds. 

Figure  2  shows  the  depth  error  verses  origin  time  error,  with  respect  to  those  of 
ISC.  The  new  points  should  show  a  linear  trend  based  on  the  assumed  earth  model. 
The  least  square  fit  of  the  data  yields: 

T=0.141jy-1.43  (6) 

where  T,  origin  time  error  in  seconds,  H  depth  mislocadon  in  km  is  equal  to  relocated 
source  depth  minus  ISC  source  depth.  The  standard  error  is  0.55  seconds  if  we 
assume  that  the  origin  times  do  not  have  errors  and  the  standard  error  is  0.91  seconds 
if  the  error  bars  are  included.  The  least  square  fitting  method  that  includes  error  bars 
assumes  that  the  probability  distribution  in  the  intervals  of  values  between  two  ends  is 
constant,  the  imn-information  probability  density.  The  coeflBcient  of  H,  0.141,  is 
dependent  on  the  input  model  TIP,  6.16  to  6.55  kra^ec  of  source  layer,  and  the  aver¬ 
age  ray  parameter. 

The  offset  of  -1.43  seconds  comes  from  the  differences  of  earth  model,  the  sta¬ 
tion  correction  and  topographic  correction  relative  to  ISC  assumption.  For  example, 
suppose  we  pick  an  event,  August  20,  1986,  21  :  23  :  53.9  (  ISC  origin  time  ).  which 
falls  near  the  line  or  on  the  average  level  (  Figure  2  )  and  examine  the  various  correc¬ 
tions.  We  relocate  it  and  find  a  new  origin  time  21  :  23  :  54.21.  In  this  calculation,  we 
assume  the  source  depth  of  25  km  given  by  ISC,  and  the  JB  model  for  both  source 
region  and  the  receiver  structures,  without  correcting  for  the  station  delay.  The  travel 
time  difference  resulting  from  replacing  the  JB  model  by  TIP  is  0.4  seconds  assuming 
a  ray  parameter  of  0.07  sec^.  Station  correction  contributes  about  -0.2  seconds  to 


\ 
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the  origin  times.  'Die  elevation  correction  of  5  km  results  in  0.8  seconds.  The  total  of 
all  these  corrections  accounts  for  the  1.4  second  offset 

Discussion 

Barazangi  and  Ni  (  1682  ),  and  Mi  and  Barazangi  (  1683  )  used  P„-  and  5„- 
waves  crossing  the  Tibetan  Plateau,  and  obtained  velocities  of  8.42  kmyfeec  for 
and  4.73  kmyfeec  for  5„.  Holt  and  Wallace  (  1660  )  inverted  the  waveform  daUi 
and  concluded  that  the  P,  velocity  is  8.24  km/feec  beneath  the  Tibet  region.  Zhao  et  al. 
(  1661  )  reported  that  the  P„  velocity  beneath  the  Tibetan  Plateau  is  8.3  km^^ec,  and 
that  the  S,  velocity  is  4.6  km^ec  by  modeling  SS  and  S  waveforms.  Zhao  et  aJ.  ( 
1661  )  (Ed  not  use  the  absolute  travel  time  in  developing  the  velocity  model  TIP,  but 
rather  the  differential  travel  times  of  SSS,  thus  their  results  on  the  upper  mantle  velo¬ 
cities  are  not  effected  by  the  accuracy  of  the  origin  times.  Using  P^i  waveforms  and 
the  travel  times  of  P„  and  S„  phase  may  introduce  two  percent  higher  lid  velocity  bias 
due  to  the  non-planar  characteristics  of  the  Moho  boundary  beneath  the  Tibetan  Pla¬ 
teau  as  given  by  Zhao  and  Helmberger  (  1891  ). 

Models  based  on  absolute  travel  times  obviously  do  depend  on  the  origin  time. 
For  example.  Holt  and  Wallace  (  1990  )  used  24  events  from  population  given  in 
Table  1.  They  assumed  these  events  are  shallow  with  depths  similar  to  the  relocated 
depths  given  in  Table  1,  and  applied  the  ISC  origin  times  with  corrections  of  the 
source  depth,  cEpping  Moho  and  70  km  crust  If  they  did  not  correct  the  origin  times, 
their  P„  travel  times  would  be  3.6  seconds  too  short  relative  to  the  results  presented  in 
Table  1  on  average.  Although  they  over  corrected  the  effects  of  the  velocity  structure 
of  source  region,  two  seconds,  (  they  assume  that  below  the  crust  the  Tibetan  Plateau 
has  a  JB  velocity  cEstribution  ),  their  result  of  P„  velocity  8.24  km/feec  is  reasonable. 

The  origin  time  effects  on  the  P„  and  results  from  the  regional  phases  are 
slightly  less  dramatic,  if  both  the  source  depths  and  origin  times  of  ISC  are  used.  This 


87 


is  because  the  deeper  source  depth  of  ISC  compensates  its  later  origin  time  The 
effects  of  origin  time  errors  and  depth  errors  on  the  travel  times  of  P„  and  S\  are 
given  by 

Arp^=-0.034i/-i-1.43  (7) 

and 

Arj.»0.044^+1.43  (8) 

assuming  that  the  compressionaJ  velocity  6.2  km/feec,  the  shear  velocity  3.5  km  /feec 
for  the  crust  and  velocities  for  the  mantle  are  8.3  and  4.6  Arms"*  respectively.  Averag¬ 
ing  over  all  the  relocated  earthquakes,  yields  A7]b„  =  -1-1.8  seconds  and  A75„  =  -1-0.95 
seconds. 

The  travel  time  of  phase  recorded  at  a  station  10  degrees  away  from  the 
source  is  about  140  seconds.  From  Ekjuation  (7),  the  P„  travel  time  is  1.8  seconds 
longer  than  that  predicted  by  using  ISC  origin  time  and  source  depth,  averaging  over 
all  the  relocated  earthquakes.  Considering  this  travel  time  difference,  we  get  a  P„  velo¬ 
city  of  8.43  Arms"*  by  using  the  P„  travel  time  predicted  by  using  ISC  source  deptli 
and  origin  time  if  the  real  P„  velocity  is  8.30  Attis"*.  This  velocity  difference  almost 
covers  the  possible  models  proposed  for  the  Tibet  region  discussed  earlier. 

The  travel  time  of  a  source-receiver  cDstance  of  10  degrees  is  about  260 
seconds.  From  Equation  (8),  the  5„  travel  time  is  0.95  seconds  longer  than  that 
predicted  by  using  ISC  origin  time  and  source  depth,  averaging  over  all  the  earth¬ 
quakes  we  relocated.  Considering  this  travel  time  difference,  we  get  a  velocity  of 
4.62  Arms"*  by  using  5„  travel  time  predicted  by  using  ISC  source  depth  and  origin 
time  if  the  real  5,  velocity  is  4.60  Aros"*.  Thus  the  origin  time  and  source  depth 
given  by  ISC  do  not  make  very  much  difference  to  the  velocity  if  only  S„  travel 
times  are  used.  In  short  the  travel  time  errors  introduced  by  the  ISC  origin  times  and 
source  depths  lead  to  about  1.5  per  cent  over-estimate  of  the  P„  velocity  and  about  0.5 
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per  cent  over-estimated  of  the  S„  velocity. 

The  ISC  location  of  the  earthquakes  are  not  badly  mlslocated.  However  this 
mislocation  may  introduce  more  than  5  seconds  to  regional  Love  wave  travel  times, 
since  the  group  velocity  of  the  maximum  amplitude  of  Love  waves  is  about  3.5 
km/hec.  Thus,  this  small  mislocation  can  lead  to  a  bias  in  the  velocity  model.  For 
example,  suppose  we  consider  the  events  discussed  by  Zhao  et  al.  (  1991  ),  see  Table 
2.  In  this  table,  ^  is  the  distance  between  the  old  and  new  location  in  kms  at  a 
number  of  stations  surrounding  Tibet  Positive  values  mean  that  the  new  distances  are 
larger,  is  the  synthetic  travel  time  correction  caused  by  the  distance  change 

assuming  a  3.5  km/feec  group  velocity.  is  the  observed  travel  time  change  due  to 

the  change  in  origin  time,  and  6T  is  defined  as  T'obt~'^»yn‘  From  this  table,  we  see  that 
the  velocity  structure  derived  by  using  the  travel  times  of  Love  waves  are  about  one 
percent  exaggerated.  Zhao  et  al.  (  1991  )  claim  that  TIP’s  crustal  model  is  mainly  for 
the  southern  part  of  the  Tibetan  Plateau,  If  we  include  the  effects  of  the  new  origin 
times  and  locations  given  here,  the  crustal  velocity  slnicture  appropriate  for  the  north¬ 
ern  part  of  the  Rateau  is  about  three  percent  slower  than  Indicated  by  TIP. 

In  conclusion,  we  investigated  more  than  90  per  cent  of  the  Tibetan  earth¬ 
quakes  and  found  that  most  of  these  events  are  shallower  than  reported  by  ISC;  most 
are  no  deeper  than  20  km.  We  relocated  59  of  them,  and  found  that  the  average  origin 
time  is  3.06  seconds  earlier  than  that  reported  by  ISC.  Eleven  of  them  are  more  than 
5.5  seconds  earlier  than  reported.  Applying  these  corrections  to  and  travel  times 
explains  why  many  authors  have  over-estimated  the  upper  mantle  velocities  beneath 
the  Tibetan  Plateau. 
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Table  1.  Relocation  of  Tibetan  Earthquakes 


Date 

Origin  Time 

AT 

(sec) 

Location  (  ®N,  °E  i 

D 

Depth 

Mb 

.. 

Mo  Da  Yr 

h  m 

relc 

ISC 

relc 

ISO 

(kru) 

re>c  ISC 

03  16  64 

0105 

13  S±1  2 

198 

■nRi 

37  11  95  60 

36  95  95  50 

20  1 

10 

50 

59 

21 

09  26  64 

0046 

-31^13 

02  6 

-5.7 

29  91  80  55 

29  96  80  46 

10  6 

18 

50 

58 

32 

10  21  64 

2309 

15  8^1.3 

19  0 

-3.2 

28  22  93  87 

28  04  93  75 

106 

15 

37 

59 

32 

01  12  65 

1332 

22.1i:l.l 

241 

-20 

27.41  87.85 

27  40  87  84 

18 

15 

23 

58 

27 

02  05  66 

1512 

27.2±10 

32  9 

-57 

26  29  103.20 

26  22  103  21 

83 

4 

42 

56 

39 

02  13  66 

1044 

36  0±12 

380 

-20 

26.27  103  25 

26  17  103  25 

10  6 

4 

6 

54 

32 

03  06  66 

0210 

52  1±0  8 

520 

+0.1 

31  50  80  53 

31  51  80  55 

22 

11 

5 

53 

33 

03  06  66 

0215 

49  Si;  11 

57.2 

-74 

31.48  80  50 

31  49  80  50 

8 

50 

60 

37 

06  27  66 

1041 

04  2i;l  2 

08  1 

-39 

29  70  80  89 

29  62  80  83 

15 

33 

60 

42 

08  15  66 

0215 

29  9±1  1 

280 

+19 

28  70  78  91 

28  67  78  93 

45 

25 

5 

56 

31 

09  26  66 

0510 
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Table  2. 
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Figure  1.  The  source  depths  before  and  after  relocatjon.  The  top  plot  shows  the  loca¬ 
tions  of  the  earthquakes  whose  depths  are  plotted  below.  The  dote  are  the  earth¬ 
quakes  on  the  Tibetan  boundary,  the  stars  are  Inside  Tibet  The  middle  plot  shows 
the  source  depth  distribution  after  relocation  and  the  bottom  plot  shows  ISC 
depths.  The  solid  symbols  indicate  the  events  that,  we  relocated  in  this  study. 


Origin  Time  Error  (  sec 


Depth  Error  (  km  ) 


Figure  2.  Origin  time  error  verses  depth  error  respects  to  the  ISC  reports.  Ihe  vertical 
fines  show  the  error  bars. 
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